S  AD-A092  149  NAVAL  RESEARCH  LAB  WASHINGTON  DC  F/G  20/11 

EFFECTIVE  SHEAR  MODULUS  FOR  FLEXURAL  AND  EXTENSIONAL  WAVES  IN  A— ETCCU) 
SEP  80  P  S  DUBBELOAY 

UNCLASSIFIED  M«L-fl372 _ NL _ _ 


NRL  Report  8372 


Effective  Shear  Modulus  for  Flexural  and  Extensional 
Waves  in  an  Unloaded  Thick  Plate 


Pieter  S.  Dubbelday 

Underwater  Sound  Reference  Detachment 
P.  0.  Box  8337 
Orlando,  Florida  32856 


September  29,  1980 


& 


NAVAL  RESEARCH  LABORATORY 
Washington,  D.C. 

Approved  for  public  release;  distribution  unlimited. 


DTIC 

^ELECTEES 

V-  * 

\  MOV  2  6  19S0  U 


80  11  25  035 


ii’iikt* 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (#*h*n  Data  Enlarad) _ 

REPORT  DOCUMENTATION  PAGE  before^completSg^orm 

)l  REPORT  NUWB'TH"'- - ;  |2.  GOVT  ACCESSION  NO.  3  RECIPIENT>-«ATALOO  NUMBER 


Cl  REPORT  N L 


NRL  WSpwt-8372 


)-fi09A 


4  TITLE  fad  Subthit)'  _  .  5-  TYPE  OF  REPORT  ft  PERIOD  COVERED 

k‘ '  . .  . .  *7  Interim  report  on  a 

£FFECTIVE£HEAR>MODULUS  FOR  FLEXURAL  AND  continuing  NRL  Problem 
EXTENSION AL  WAVES  IN  AN  UNLtfSDED^THICK  / 

PLATE  •  ~  _ J 


6.  PERFORMING  ORG.  REPORT  NUMBER 


Pieter  S.^Dubbelday 


"J  77‘i 


/•  «£.. 


19.  PERFORMING  ORGANIZATION  NAME  AND  AODRESS 


i  •  CONTRACT  OR  GRANT  NUMBER^) 


N00173-D-0007 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  *  WORK  UNIT  NUMBERS 


61153N;  RRA110842; 
82-0585-0-0 


Underwater  Sound  Reference  Detachment  , —  61153N;  RRA1 10842; 

Naval  Research  Laboratory  //*?'<\  Zr  J  I  82-0585-0-0 

P.O.  Box  8337,  Orlando,  FL  32856  /  /oC  )  p  fa  / _ 

n.  CONTROLLING  OFFICE  NAME  ANO  AODRESS  ^ - -  - — — *  12-  REPORT  DATE 

^  ^  /  '  /  •  September  29, 1980 

^  /6  /Y  h  (7  -2j/  ^NUMBER  OF  PAGES 

I^MorflTORlNG  AGENCY  NAME  A  ADDRESS(*/  dlltarant  (root  Controlling  Olflca)  15.  SECURITY  CLASS,  (o t  thia  report) 


III.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


/6  ‘  A  t 


V  /  "  / 


vT  /  .*  /-'-v 


UNCLASSIFIED 


15a.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


|16-  DISTRIBUTION  STATEMENT  (of  thia  Raport) 


Approved  for  public  release;  distribution  unlimited. 


I  17.  DISTRIBUTION  STATEMENT  (of  tha  abstract  an  tar ad  In  Block  20,  If  dlf/arant  from  Raport) 


tS-  SUPPLEMENTARY  NOTES 

♦The  author  is  Professor  of  Physics  and  Oceanography  at  the  Florida  Institute  of 
Technology  in  Melbourne,  Florida,  working  with  Oak  Harbour  Marine  Assoc.,  ttcon- 
tractor  of  the  Underwater  Sound  Reference  Detachment  of  the  Naval  Research  Laboratory. 

19.  KEY  WORDS  (Contlnua  on  ravaraa  atda  If  nacaaaary  and  Idantlfy  by  block  numbar) 


20. '^BSTRACT  (Confirm#  on  ravaraa  atda  It  nacaaaary  and  Idantlfy  by  block  numbar ) 

Thin-plate  theory  suffers  from  failure  to  properly  predict  the  correct  phase  speed  of 
bending  waves  and  extensionaJ  waves  in  the  limit  of  high  frequency  or  large  thickness. 
Mindlin  has  adapted  the  basic  ideas  that  Timoshenko  developed  for  bars  to  the  case  of 
bending  waves  in  plates  to  partly  correct  this  situation.  Kane  and  Mindlin  have  given  a 
similar  correction  for  extensions!  waves  in  plates.  In  this  report  a  new  derivation  of  such 
correction  factors  is  presented  that  is  based  on  a  comparison  of  the  approximate  theory  , 

(Continued)  \_; 


DD  i^n,  1473  EDITION  OP  I  NOV  «5  IS  OBSOLETE 

S  01 02- L  F -01 4-6601 


ILASSihit1 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (•BNfb  Dai*  BnlAfpF; 


-  // 


'Y//- 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (WH,n  Dmlm  Enffd) 


Abstract  (Continued) 

with  the  results  of  exact  elasticity  theory  (Lamb  waves).  Explicit  equations  are  given 
to  compute  the  effective  shear  modulus  in  antisymmetric  waves.  It  is  shown  analytically 
that  the  phase  speed  calculated  with  these  correction  factors  for  the  shear  modulus 
asymptotically  approaches  the  Rayleigh  wave  speed  for  high  frequency.  The  same 
principle  of  comparison  of  approximate  theory  with  exact  elasticity  theory  is  applied 
to  the  other  terms  in  the  equations  of  motion  besides  those  depending  on  the  shear 
modulus.  This  has  not  been  entirely  successful,  but  it  promises,  among  other  things, 
to  improve  the  identification  of  the  second  root  in  the  quadratic  dispersion  relation 
with  the  first  order  antisymmetric  and  symmetric  Lamb  waves. 


Accession  For  ~  "j 

NT1S  GR.A&I  *  1 

DTIC  TAR  • 

Unannounced  j 
o  us  t  .i  /.  i  c ci  t  i  on.  _ 


i  •  y  / 

j  Avail;..'  ;  n  w  Codes 

I  ■■  i.  and/ or 

!  D  i  a.  t.  I  L^oinl 


DTK 


U. 1  ■  i 


\trv  NOV  n  o  »3C 


SECUHiTV  CLASSIFICATION  DERMIS  P AGEf»**n  Df 


CONTENTS 

List  of  Symbols . v 

INTRODUCTION . 1 

THEORY  OF  ELASTICITY  APPLIED  TO  WAVE 
PROPAGATION  IN  PLATES . 3 

Theory  of  Elasticity . 3 

Lamb  Waves . 6 

APPROXIMATE  THEORIES  FOR  WAVE  PROPAGATION  IN  PLATES  ....  9 

Thin-Plate  Theory:  Bending  (Flexural)  Waves  .  .  .  - . 10 

Thin-Plate  Theory:  Extensional  Waves . 12 

PLATE  STRESS  EQUATIONS  OF  MOTION . 13 

Antisymmetric  Waves . 13 

Symmetric  Waves . 14 

THICK-PLATE  THEORY . 15 

Timoshenko-Mindlin  Plate  Theory . 16 

Symmetric  Waves  in  Thick  Plates . 20 

EFFECTIVE  SHEAR  MODULUS  IN  ANTISYMMETRIC  WAVES . 23 

EFFECTIVE  SHEAR  MODULUS  IN  SYMMETRIC  WAVES . 26 

CONCLUSIONS  AND  PLANS  FOR  FURTHER  WORK . 30 

ACKNOWLEDGMENTS . 31 

REFERENCES . 31 

APPENDIX  A  —  Formulae  for  Field  Variables  from 

Elasticity  Theory . 33 

APPENDIX  B  —  Generalized  Correction  Factors . 37 

APPENDIX  C  —  Quadratic  Approximation  of  Displacement  Components  ...  40 


APPENDIX  D  —  Comparison  With  Other  Work  on  Approximate  Plate  Theory  .  .  43 


LIST  OF  SYMBOLS 


Ag  amplitude  of  symmetric  part  of  potential  <j> 

Ba  amplitude  of  antisymmetric  part  of  potential  <t> 

Ca  amplitude  of  antisymmetric  part  of  potential  \p 

c  phase  speed 


cb 

cd 


phase  speed  of  bending  waves;  c2b  = 
phase  speed  of  dilational  waves; 


E(kd)2 

3/t;(l-v2) 

E(  1-v) 

'  p(l+^)(l-2^) 


c  phase  speed  of  extensional  waves;  c2  _  E 

P(l- 1'2) 

cR  phase  speed  of  Rayleigh  waves 

E 

c „  phase  speed  of  shear  waves;  cl  =  - - 

s  s  2p(l+^) 

d  one-half  of  plate  thickness 

2  Fd 2 

D  bending  stiffness  of  plate;  D  =  - 5— 

3(l-r2 ) 

Dg  amplitude  of  symmetric  part  of  potential 

e  =  e  +  ev  +  e,;  dilation 

E  Young’s  modulus 

f  frequency 

G  shear  modulus  (=p) 

G'  effective  shear  modulus  in  antisymmetric  waves 

G"  effective  shear  modulus  in  symmetric  waves 


v 


r 


k  wavenumber 

kd  wavenumber  of  dilational  wave 

ka  wavenumber  of  shear  wave 

Lx  =  fzudz;  moment  of  u  displacement 

Ly  =  fzvdz ;  moment  of  v  displacement 

Lz  =  fzwdz;  moment  of  w  displacement 

Mx  =  fzoxdz;  bending  moment 

Af  =  fzorvdz ;  twisting  moment 

My  =  fzOydz  -,  bending  moment 

N  =  foxdz\  stress  integral 

Nv  =  Jo  dz\  shear  stress  integral 

Ny  =  foydz;  stress  integral 

Nz  =  fozdz;  stress  integral 

q  =  ( k2-k ^ 

Qx  =  f°zxdz\  shear  stress  integral 

Qy  =  fozydz\  shear  stress  integral 

r  coefficient  of  quadratic  form  in  expansion  of 

Rx  =  fzazxdz\  shear  stress  moment 

R  =  fzo  dz;  shear  stress  moment 

y  zy 

s  ( k2-k])l/i 

s  displacement  vector 

sd  irrotational  part  of  s 

$s  solenoidal  part  of  s 

Tx  =  fudz;  u  displacement  integral 

Ty  =  fvdz;  v  displacement  integral 


vi 


=  fwdz;  w  displacement  integral 


Tz 

u,l\w  components  of  displacement  vector  s 

U,V,W  components  of  zero  order  expansion  of  s  in  terms  of  z 
l-2v 

a  (c»/c<*>  =  2(H 


y 


ib 


*>d 


7r 

e 

K? 

«3 


«4 

^6 

X 

p 

v 

P 

o 

<f> 


c/c. 


Cjf/ct 

strain  tensor 

correction  factor  for  shear  modulus  in  antisymmetric  waves 

correction  factor  for  shear  modulus  in  symmetric  waves 

correction  factor  for  bending  stiffness  and  rotational  inertia  in  antisymmetric 
waves 

correction  factor  for  stiffness  modulus  in  x-direction  in  symmetric  waves 

correction  factor  for  stiffness  modulus  in  z-direction  in  symmetric  waves 

first  Lame  constant 

second  Lame  constant  (=G) 

Poisson’s  ratio 

density  of  solid 

stress  tensor 

d<t>x  d$y 
dx  +  dy 


vii 


1 

<p  potential  for  irrotational  part  of  displacement  vector 

<t>x.4>y  components  of  rotation  angle  of  plate  cross  section 
bU  bV 

* 

♦  potential  for  solenoidal  part  of  displacement  vector 

SI  rotational  displacement  vector;  =  curl  s  =  curl  sg 

co  =  2 nf;  angular  frequency 

\p  y  component  of  ¥  for  2-dimensional  displacement 

X  coefficient  of  first-order  approximation  of  the  vertical  displacement  w  in  terms 

of  the  coordinate  z 


viii 


EFFECTIVE  SHEAR  MODULUS  FOR  FLEXURAL  AND 
EXTENSIONAL  WAVES  IN  AN  UNLOADED  THICK  PLATE 


INTRODUCTION 

The  propagation  of  straight-crested  waves  in  flat  plates  is  an  important  acoustical 
and  mechanical  problem.  The  simplest  of  such  problems  is  that  of  a  plate  in  vacuum, 
where  the  plate  thickness  is  small  compared  with  the  wavelength  of  a  wave  propagating 
parallel  to  the  faces  of  the  plate.  In  this  case  the  phase  speed  of  two  types  of  waves  can 
be  readily  derived,  namely  bending  or  flexural  waves  and  extensional  or  quasi-longitudinal 
waves.  The  theory  explaining  these  waves  is  indicated  by  thin-plate  theory  or  classical- 
plate  theory.  On  the  other  hand  one  also  finds  application  of  exact  theory  of  elasticity 
to  plates  without  restriction  to  small  thickness,  which  leads  to  the  so-called  Lamb  waves. 

The  simplicity  of  thin-plate  theory  as  compared  with  the  greater  complexity  of  exact 
theory  prompts  one  to  look  for  corrections  in  those  cases  where  thin-plate  theory  gives 
unacceptable  results.  Specifically,  at  high  frequencies  (or  greater  thicknesses),  the  phase 
speed  of  bending  waves  grows  beyond  bounds  and  the  phase  speed  of  extensional  waves 
stays  constant.  Thus  neither  speed  approaches  the  value  for  Rayleigh  surface  waves,  as 
one  might  expect  on  physical  grounds. 

To  remedy  this  situation,  Mindlin  [1]  adapted  the  ideas  of  Timoshenko  for  propa¬ 
gation  along  bars  to  the  case  of  plates,  namely  the  introduction  of  the  effects  of  trans¬ 
verse  shear  stress  and  rotatory  inertia.  In  his  development,  Mindlin  introduces  an  effec¬ 
tive  shear  modulus  to  account  for  the  difference  between  the  simplified  shear  angle  profile 
profile  and  the  actual  one.  The  ensuing  corrective  factor  is  fixed  in  such  a  way  that  the 
phase  speed  approaches  the  Rayleigh  wave  speed  asymptotically  for  high  frequencies. 

Although  Mindlin  succeeds  by  this  manner  in  devising  a  viable  alternative  to  the 
Lamb  wave  solution,  without  the  drawbacks  of  thin-plate  theory,  there  is  a  certain  un¬ 
satisfactory  aspect  to  the  determination  of  the  correction  factor  in  the  sense  that  its 
value  is  imposed  from  outside  the  theory  proper.  This  has  been  mentioned  by  other 
authors,  for  example  Cremer  et  al.  [2].  Cowper  [3]  also  discusses  various  criticisms  and 
solutions  for  the  similar  issue  in  the  case  of  bars.  In  applications  of  the  Mindlin-Timo- 
shenko  plate  theory,  it  is  not  always  appreciated  that  the  introduction  and  evaluation 
of  a  correction  factor  depend  on  the  physical  circumstances  of  the  problem,  and  the 
factor  is  certainly  not  a  universal  constant.  As  an  example,  in  a  paper  on  wave  propa¬ 
gation  in  a  plate  loaded  with  an  incompressible  fluid  layer  (surface  waves  only),  Walter 
and  Anderson  [4]  state  that  the  pertinent  correction  factor  is  the  root  of  the  dispersion 
relation  for  Rayleigh  waves  for  a  solid  in  vacuum.  However,  the  dispersion  relation  for  a 
fluid-loaded  plate  is  not  the  same  as  for  a  plate  in  vacuum. 

This  state  of  affairs  prompted  the  question  of  whether  it  would  be  feasible  to  de¬ 
sign  a  new  method  of  determining  the  correction  factor  based  on  the  way  it  is  intro¬ 
duced  into  the  theory.  Instead  of  choosing  the  factor  in  such  a  way  that  it  reproduces  the 
Rayleigh  wave  speed  at  high  frequency,  as  was  done  by  Mindlin  [1],  this  method  should 
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display  the  proper  frequency  dependence  of  the  factor  in  such  a  way  that  at  high  fre¬ 
quencies  the  Rayleigh  wave  speed  emerges  automatically.  Once  a  method  of  computation 
is  found  for  a  correction  factor  operating  properly  for  the  simple  case  of  a  plate  in 
vacuum,  it  is  expected  that  this  same  method  can  be  applied  to  more  complicated  cases: 
plates  loaded  by  a  fluid,  on  one  or  both  sides;  composite  plates,  possibly  with  fluid 
loading. 

Originally  the  emphasis  in  the  analysis  was  placed  on  the  case  of  antisymmetric 
waves  only,  in  connection  with  Timoshenko-Mindlin  theory.  It  appears  that  these  waves 
are  more  often  referred  to  in  hydroacoustical  applications  than  symmetric  waves.  When 
a  sound  wave  in  a  fluid  impinges  on  a  plate,  both  symmetric  and  antisymmetric  waves 
are  induced  and  it  is  necessary  to  know  how  the  acoustical  energy  is  divided  over  the 
two  wave  types  in  order  to  calculate  the  acoustic  radiation,  reflection,  or  absorption. 
Therefore,  while  embarking  on  a  new  evaluation  procedure  for  the  correction  factor  in 
the  Timoshenko-Mindlin  theory,  it  became  clear  that  a  similar  procedure  should  be 
applied  to  the  extension  of  the  theory  for  extensional  waves  to  greater  plate  thichnesses 
(or  higher  frequencies).  An  analysis  of  this  type,  from  a  different  viewpoint,  however,  was 
performed  by  Kane  and  Mindlin  [5].  Their  discussion  does  not  stress  the  basic  analogy 
between  the  antisymmetric  and  symmetric  waves  that  is  an  important  feature  of  the 
present  study.  A  previous  example  of  considering  Timoshenko-Mindlin  and  Kane-Mindlin 
theories  simultaneously  is  found  in  the  article  by  Walter  and  Anderson  [4],  referred  to 
above.  Their  study  is  of  limited  direct  value  for  hydroacoustics,  since  the  fluid  loading  the 
plate  is  assumed  incompressible  and,  as  a  consequence,  only  surface  waves  are  admitted. 
These  authors  do  not  give  an  independent  evaluation  of  the  correction  factors. 

In  the  present  study,  the  first  section  gives  a  short  tutorial  introduction  to  the  sub¬ 
ject  of  elasticity  theory  and  Lamb  waves.  This  is  followed  by  a  discussion  of  thin-plate 
theory,  with  the  explanation  of  flexural  and  extensional  waves.  The  subsequent  deriva¬ 
tion  of  plate  stress  equations  of  motion  is  still  mainly  tutorial,  but  emphasizes  the  analogy 
between  antisymmetric  and  symmetric  waves  and  their  differences,  based  on  the  parity 
of  the  field  variables.  This  approach  leads  to  a  rationale  for  the  two  types  of  waves.  The 
section  on  thick-plate  theory  describes  the  work  of  Mindlin  for  antisymmetric  waves  and 
a  parallel  treatment  for  the  case  of  symmetric  waves. 

The  comparison  of  the  thick-plate  equations  of  motion  with  the  corresponding  equa¬ 
tions  in  exact  elasticity  theory  leads  to  a  method  of  calculating  a  correction  factor  for 
shear  modulus  both  in  antisymmetric  and  symmetric  waves.  The  resulting  phase  speeds 
possess  the  property  desired  from  a  viewpoint  of  physical  intuition,  namely  that  they 
asymptotically  approach  the  phase  speed  of  Rayleigh  surface  waves  at  high  frequencies. 
Detailed  formulae  for  the  field  variables  from  elasticity  theory  are  given  in  Appendix  A. 

It  would  appear  logical  to  extend  this  same  approach  to  a  comparison  of  other  ap¬ 
proximate  terms  in  the  equations  of  thick-plate  theory  with  those  of  elasticity  theory 
and  to  establish  additional  correction  factors.  Certain  difficulties  were  encountered  in 
this  development  and  because  of  the  tentative  character  of  this  theory,  it  is  discussed 
in  appendices  B  and  C  to  the  main  body  of  the  report.  These  other  correction  factors 
are  helpful  in  the  identification  of  the  second  root  in  the  quadratic  dispersion  relation 
with  the  first  mode  in  the  antisymmetric  and  symmetric  Lamb  waves. 
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In  the  vast  literature  on  plate  theory  there  are  various  articles  and  reports  that 
show  an  outward  similarity  to  the  work  reported  here.  The  distinction  between  those 
papers  and  the  present  analysis  may  be  difficult  to  appreciate.  Therefore  in  Appendix 
D,  a  discussion  is  presented  of  the  relation  of  the  work  reported  here  to  other  work 
in  approximate  plate  theory. 


THEORY  OF  ELASTICITY  APPLIED  TO  WAVE  PROPAGATION  IN  PLATES 
Theory  of  Elasticity 

The  basic  equations  of  elasticity  theory  are  derived  in  various  texts  (see  e.g.  Timo¬ 
shenko  and  Goodier  [6]).  The  relationships  between  the  elements  of  the  stress  tensor  a 
and  the  strain  tensor  e  are  given  by 


where 


ax=  \e  +  2  Gex 

°xy 

Ge  xy 

ay  =  \e  +  2  Gey 

°yz 

=  Geyz 

(1) 

az  =  Xe  +  2  Gez 

°zx 

=  Gezx, 

X  =  first  Lame  constant 

G  =  shear  modulus  (=  p  second  Lame  constant) 
e  =  dilation  =  ex  +  ey  +  ez. 

The  elements  of  the  strain  tensor  are  given  in  terms  of  the  displacement  vector  s  of  a 
particle  with  components  u,  v,  w  by 


and 


du 

du 

dv 

x  dx 

€xy 

dy 

+  dx 

dv 

dv 

dw 

y  dy 

ey* 

dz 

+  ay 

d  w 

dw 

du 

1  N 
I'-o 

1 

N 

€zx 

dx 

dz 

(2) 
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The  equations  of  motion  in  terms  of  ihe  stress  components  are 


3o 

dor„ 

dci 

d  2  u 

X 

+ _ *y 

+  xz  -  p 

dx 

dy 

bz 

dt2 

d(JyZ 

do 

+  *  =  Q 

d2u 

dx 

3y 

dz 

dt2 

dozx 

+  3a*y 

do. 

+  _ _  =  n 

b2  w 

dx 

3y 

dz 

dt2 

where  p  is  the  density  of  the  solid. 


(3) 


Like  any  vector  field,  the  displacement  field  s  can  be  represented  as  the  sum  of  an 
irrotational  field  sd  and  a  solenoidal  field  ss  (Helmholtz  representation).  Therefore, 

s  =  sd  +  sg,  (4) 

where  curl  sd  =  0  and  div  sg  =  0.  As  a  consequence,  for  the  dilation, 

e  =  div  s  =  div  sd;  (5) 

and  for  the  rotational  displacement, 

£1  =  curl  s  =  curl  sg.  (6) 

Waves  related  to  the  variations  of  sd  are  indicated  by  dilatational  waves;  those  con¬ 
nected  with  sg  are  called  shear  waves.  Only  in  media  of  infinite  extent  do  these  waves 
occur  in  a  pure  form:  in  finite  media  the  two  types  mix  as  a  consequence  of  the  boundary 
conditions.  In  fluids  only  dilatational  waves  can  occur.  The  terms  longitudinal  and  trans¬ 
verse,  respectively,  have  to  be  used  with  caution;  a  dilatational  wave  is  longitudinal  only 
if  it  is  a  plane  wave,  and  a  shear  wave  is  transverse  only  if  it  is  a  plane  wave. 


The  displacements  in  a  dilatational  wave  can  be  derived  from  a  scalar  potential  <p, 
where 


sd  =  grad  0,  (7) 

and  the  displacements  in  a  shear  wave  can  be  derived  from  a  vector  potential  if* ,  where 

sg  =  curl  0.  (8) 

Whenever  one  has  a  2-dimensional  wave  propagation,  only  one  component  of  the 
potential  0  is  nonzero.  This  component  is  indicated  by  0. 
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For  the  case  of  a  dilatational  wave  the  dilatation  e,  the  displacement  vector  s,  and 
the  potential  t!/  are  all  solutions  to  the  differential  equation 

a  2 

(X+2G)  v  2  \e,s4  \  =  9  (9) 

Similarly  for  a  shear  wave,  the  relevant  wave  equation  for  the  rotational  displacement 
fl,  the  displacement  s,  and  the  vector  potential  t/i  is 

d2 

g  2  \  n,s,+\  =  P—\n,s,t\.  (10) 

Notice  in  the  above  that  the  symbolic  operator  2,  when  operating  on  a  vector,  is  given 
only  by  the  Laplacian  form 


(11) 


for  Cartesian  coordinates.  For  curvilinear  coordinates  this  operator  should  be  interpreted 
as 


v2  =  grad  div  -  curl  curl.  (12) 

The  wave  equations  show  that  dilatational  waves  propagate  with  speed 

cd  =  [(X+2G)/p):/2,  (13) 

where  cd  is  the  phase  speed  of  dilatational  waves.  Also,  shear  waves  propagate  with  speed 

cs  =  (G/p)V\  (14) 

where  cs  is  the  phase  speed  of  shear  waves. 

It  may  be  pointed  out  here  that  in  general  the  elastic  properties  of  solids  may  be 
represented  by  two  independent  constants  for  which,  thus  far,  X  and  G  have  been  used. 
Another  useful  set  is  Young’s  modulus  E  and  Poisson’s  ratio  v.  The  following  conversion 
will  be  extensively  used  below.  The  shear  modulus  is  given  by 

<i5i 


and  the  “stiffness  modulus”  by 


X  +  20  = 

(1+p)(1-2p)  ' 


(16) 
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Lamb  Waves 

Lamb  waves  refers  to  the  2-dimensional  straight-crested  plane  waves  in  a  homogeneous, 
isotropic  plate  in  the  direction  parallel  to  the  faces  of  the  plate  (see  Fig.  1  for  the 
geometry).  Notice  that  the  thickness  of  the  plate  equals  2d.  Some  authors  use  the  symbol 
h  for  the  plate  thickness,  e.g.,  Mindlin  [1J. 


Z 


Fig.  1  —  Plate  geometry 


The  derivation  of  the  theory  of  Lamb  waves  given  below  follows  clo*  :he  develop¬ 
ment  of  Viktorov  [  7a] . 

The  waves  can  be  most  advantageously  expressed  in  terms  of  the  potentials  0  and  0 
in  the  form 


0  =  t^4s  cosh(qz)  +  Ba  sinh(qz)]  exp[i(o>f  fex)] 

(17) 

0  =  [Ds  sinh(sz)  +  Ca  cosh(sz)]  exp[i(cot  fex)] 

The  wave  propagation  velocity  c,  the  phase  speed,  equals  w/fc,  where  u>  is  the  angu¬ 
lar  frequency,  k  is  the  wavenumber,  and  the  symbols  q  and  s  are  given  by 

q 2  =  k2-k2d  =  k2  ji-(c/cd)2| 

and  (18) 


s2  =  k2-k\  =  k2  |l-(c/cg)2|, 

where  kd  and  k$  are  the  dilatational  and  shear  wave  numbers,  respectively, 

kd  =  cj[p/(X+2G)]'/2  =  k(c/cd) 


and 


(19) 


ks  =  w(p/G)1/4  =  k(c/cs)  . 

The  subscripts  a  and  s  in  the  amplitudes  46,  Ba,  Ds,  and  C0  refer  to  the  two  possible 
types  of  Lamb  waves,  antisymmetric  and  symmetric.  These  terms  indicate  the  symmetry 
character  of  the  cross  section  of  the  plate  (see  Fig.  2).  Parity  of  the  displacement  func¬ 
tions  u  and  w  with  respect  to  the  z  coordinate  is  different  for  the  two  types:  in  anti¬ 
symmetric  waves  u  is  odd  and  w  is  even,  whereas  in  symmetric  waves  u  is  even  and  w 
is  odd.  The  displacement  components  are  derived  from  the  potentials  by 
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30  30 

dx  3  z 


and  (20) 

3  0  30 

w  =  3T  "  37' 


This  explains  the  various  products  of  the  amplitudes  and  hyperbolic  functions  in  Eq.  (17). 


ANTISYMMETRIC 


SYMMETRIC 


Fig.  2  —  Lamb  waves 

The  amplitudes  are  related  through  the  boundary  conditions,  which  are  the  values 
for  the  normal  and  shear  stresses  applied  at  the  two  faces  of  the  plate.  These  stresses  are 
found  in  terms  of  the  amplitudes  by  means  of  Eqs.  (1)  and  (2).  This  leads  to  the  follow¬ 
ing  set  of  expressions  for  the  stresses,  using  the  definitions  of  Eq.  (18); 

az  ~  G|>ls(ft2+s2)cosh(qz)  +  Bfl(k2+s2)sinh(gz) 

-  Ca  2 iks  sinh(s2)  -  Ds  2iks  cosh  (,sz)J 

(21) 


azx  =  -G|/ls  2 ikq  sinh(<?2)  +  Ba  2 ikq  cosh(<?2) 

+  Ca(fc2+s2)cosh(sz)  +  Z?s(fc2+s2)sinh(sz)|. 

The  applied  stresses  for  a  plate  in  vacuum  are  zero.  In  that  case  the  antisymmetric  and 
symmetric  waves  independently  satisfy  the  zero  boundary  conditions,  and  thus 

J30(fe2+s2)sinh(q|d)  -  C0  2ifes  sinh(sd)  =  0 

Ba  2 ikq  cosh (qd)  +  Ca(k2+s2)cosh(sd)  =  0 

Ag(k2+s2)cosh(qd)  -  Ds  2 iks  cosh (sd)  =  0 


and 


(22) 


As  2ikq  sinh(gd)  +  Dg(fe2+s2)sinh(sd)  =  0. 
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By  elimination  of  the  amplitudes  from  these  equations,  one  arrives  at  the  dispersion 
relations 

(fe2+s2)2  sinh(qd)cosh(sd)  -  4 k2qs  cosh(qd)sinh(sd)  =  0 
for  antisymmetric  waves  and 

(fe2+s2  )2cosh(gd)sinh(sd)  -  4k2 qs  &mh(qd)cosh(sd )  =  0  (24) 

for  symmetric  waves.  At  large  values  of  fed  both  dispersion  relationships  approach  the 
same  relation, 

4k2 qs  -  (fe2+s2)2  =  0,  (25) 

which  is  the  dispersion  relation  for  Rayleigh  waves  defined  as  waves  at  the  surface  of  a 
semi-infinite  solid.  This  is  understandable  since,  for  a  plate  thickness  that  is  large  com¬ 
pared  with  the  wavelength,  the  Lamb  wave  separates  into  two  independent  surface 
waves.  The  Rayleigh  wave  corresponds  to  the  root  of  Eq.  (25)  for  which  the  ratio  c/cg 
lies  between  0  and  1.  A  good  approximation  for  this  root  is  the  expression  (Viktorov  (7b)) 


c  0.87  +  1.12  v 


This  ratio  varies  from  0.87  to  0.95  when  Poisson’s  ratio  v  varies  from  0  to  0.5. 

For  a  fixed  value  of  the  thickness  2d  and  the  frequency  f,  the  dispersion  relations  for 
Lamb  waves  have  a  finite  number  of  real  roots,  corresponding  to  a  finite  number  of 
modes  of  propagation  in  the  direction  of  the  plate.  There  are  infinite  purely  imaginary 
roots  that  correspond  to  waves  perpendicular  to  the  faces  of  the  plate  and  that  decay  or 
increase  exponentially  parallel  to  the  faces  of  the  plate.  Figures  3  and  4  show  the  dimen¬ 
sionless  wave  speed  c/c,  as  a  function  of  the  dimensionless  wavenumber  kgd  for  the 
zero-order  antisymmetric  and  symmetric  Lamb  waves,  respectively. 


Fig.  3  —  Dispersion  curve  for  zero-order  antisymmetric  Lamb  wave, 
with  v  =  0.355  and  7R=  0.936 
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Fig.  4  —  Dispersion  curve  for  zero-order  symmetric  Lamb  wave,  with 
v  =  0.355  and  >R=  0.936 

APPROXIMATE  THEORIES  FOR  WAVE  PROPAGATION  IN  PLATES 

Application  of  the  theory  of  elasticity  to  wave  propagation  in  plates  leads  to  rather 
complex  expressions,  expecially  when  composite  plates  and  fluid  loading  of  plates  are 
studied.  Therefore,  one  attempts  to  build  a  simpler  theory  that  makes  use  of  the  fact  that 
in  plates  one  dimension  is  much  smaller  than  the  other  two. 

For  the  case  where  the  dimensionless  wavenumber  kd  is  small,  this  theory  ex¬ 
plains  the  occurrence  of  two  types  of  wave  motions,  flexural/bending  waves  and  exten- 
sional  waves.  These  two  types  correspond  to  the  antisymmetric  waves  and  symmetric 
waves,  respectively,  in  elastic  theory.  The  phase  speed  of  bending  waves  is  proportional 
to  kd,  and  the  phase  speed  of  extensional  waves  is  constant.  Thus,  in  neither  case  does 
the  functional  dependence  on  kd  correctly  represent  the  behavior  of  the  phase  speed  at 
larger  kd,  since  it  is  plausible  on  physical  grounds  that  the  phase  speed  of  each  type  ap¬ 
proaches  the  phase  speed  of  Rayleigh  surface  waves.  To  correct  this  flaw  of  the  thin- 
plate  theory,  an  extension  of  the  theory  has  been  developed  indicated  below  by  the  term 
thick-plate  theory. 

Mathematically,  both  thin-plate  and  thick-plate  theories  can  be  comprehensively 
represented  by  the  following  series  expansion  of  the  displacement  components  in  terms 
of  the  coordinate  z  to  the  first  order  (Walter  and  Anderson  [4] ): 

u(x,y,z,t )  =  U(x,y,t)  +  z<t>x(x,y,t) 

v(x,y,z,t)  =  V(x,y,t)  +  z<fiy(x,y,t)  (27) 

ie(jc,y,a,f)  =  W(x,y,t)  +  *x(*,.v,f)- 
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The  functions  <px  and  <py  have  the  character  of  angular  rotations;  the  function  x 
represents  a  displacement  per  unit  thickness. 

Thin-plate  theory  follows  by  further  restrictions  on  the  functions  <p  ,  <t>y,  and  x  — 
restrictions  that  are  discarded  in  the  case  of  thick-plate  theory.  In  general  the  same  func¬ 
tion  symbols  will  be  used  for  time-dependent  and  time-independent  functions,  since  the 
dependence  on  time  is  assumed  to  be  only  in  exponential  form.  Thus,  for  example,  the 
time-dependent  displacement  u(x,y,z,t )  is  alternatively  expressed  as  u(x,y,z)  exp(icof). 

Thin -Plate  Theory:  Bending  (Flexural)  Waves 

Bending  or  flexural  waves  correspond  to  the  antisymmetric  Lamb  waves.  They  are 
mathematically  represented  by  the  part  of  the  expansion  Eq.  (27)  for  which  u  and  c  are 
odd  functions  of  z  and  for  which  w  is  an  even  function  of  z.  Thus, 

u(jt,y,r.{)  =  z  <px(x,y,t) 

v(x,y,z,t)  =  z  <>y(x,y,0  (28) 


and 


w(x,y,z,t)  =  W(x,v,f). 

One  might  surmise  that  bending  waves  are  predominantly  governed  by  shear  forces, 
but  the  opposite  is  ture.  the  major  elastic  force  is  due  to  compressions  and  extensions, 
as  in  the  case  for  extensions!  waves.  The  difference  with  extensional  waves  is  that  here 
the  compression  of  a  particle  above  the  neutral  plane  is  accompanied  by  an  extension 
of  a  particle  below  the  neutral  plane,  and  vice  versa.  The  thin-plate  approximation 
emphatically  excludes  shear  forces  by  making  the  assumption  that  exy  and  fyz  are 
zero,  and  thus  the  functions  <px  and  4>y  are  related  to  the  function  W  by  the  equations 

dW 


and  (29) 

dW 

~  ~  dy  ' 

Geometrically,  this  means  that  cross  sections  originally  perpendicular  to  the  neutral 
beam  will  stay  perpendicular  to  the  neutral  plane  (see  Fig.  5). 
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In  thin-plate  theory  it  is  assumed  for  both  types  of  waves  that  zero  normal  stress 
on  the  faces  of  the  plate  implies  that  o^x.y.z.t)  =  0.  Equation  (1)  shows  that  if  oz  is 
to  be  equal  to  zero,  then  one  has  (X+2G)ez  +  X(£z+ey)  equal  to  zero.  Eliminating  ez 
between  this  expression  and  the  expression  for  ax  and  a y  in  Eq.  (1)  results  in 


a 


x 


+  V  ey) 


and 


(30) 


with  the  substitutions 


4  G(\+G) 
X  +  2  G 


E  X 

'  1_„2  and  2(X+2G)  =  v 


An  essential  feature  of  approximate  plate  theory  is  that  the  equations  of  motion 
(Eq.  (3))  are  integrated  in  the  z  direction,  perpendicular  to  the  faces  of  the  plate.  De¬ 
pending  on  the  parity  of  the  displacement  functions  u,  w  with  respect  to  the  z  coor¬ 
dinate,  some  expressions  of  the  set  (Eq.  (3))  become  identically  zero  upon  integration. 
This  explains  the  different  way  in  which  the  modulus  E/(l-v2 )  of  Eq.  (30)  appears  in 
bending  waves  as  compared  with  extensional  waves.  The  stress  ox  has  odd  parity  with 
respect  to  z  for  bending  waves,  and  as  a  consequence  its  integral  over  the  plate  thick¬ 
ness  is  equal  to  zero.  Its  role  is  taken  over  by  the  moment  of  the  stress,  leading  to  the 
integral  f  +dd  zoxdz.  With  the  first  of  the  definitions  (Eq.  (2)),  the  first  expression  in 
Eq.  (30),  and  the  value  for  u  from  Eq.  (28),  this  integral  produces  a  “bending  stiffness” 

D  given  by 

2  Ed3 

D  ~  3(l-^r)  '  (3D 

Another  consequence  of  the  introduction  of  moments,  due  to  the  parity  of  the  rele¬ 
vant  functions,  is  that  the  wave  equation  for  nending  waves  is  not  of  second  order  but 
rather  fourth  order,  and  is  given  by 
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2  pd 


d2W 

at 2 


(32) 


where  the  function  W  can  be  replaced  by  any  of  the  other  field  variables.  (For  a  complete 
derivation  see,  for  instance,  Cremer  et  al.  [2b].)  By  inserting  a  traveling  wave  expression 
for  W  into  Eq.  (32),  one  finds  that  bending  waves  are  dispersive,  with  a  dispersion  relation 
for  the  wave  speed  cb, 


C2  -.E(M).y 

6  3p(l-v2 ) 


(33) 


Thin-Plate  Theory:  Extensional  Waves 

Extensional  waves  correspond  to  symmetric  Lamb  waves.  They  are  represented  by 
the  part  of  the  expansion  (Eq.  (27))  for  which  u,  v  are  even  functions  of  z  and  for  which 
w  is  an  odd  function  of  z.  Thus, 


and 


u(x,y,z,t)  =  l/(x,y,t) 

i/(x,y,2,t)  =  V(x,y,f)  (34) 

w(x,y,z,t)  =  z\(x,y,t). 


In  the  case  of  antisymmetric  waves  the  assumption  that  the  shear  stress  o  zx  is  zero  leads 
to  the  relations  of  Eq.  (29)  among  the  various  functions.  The  analogous  path  for  sym¬ 
metric  waves  leads  from  the  assumption  that  the  normal  stress  a z  is  zero  to  a  relation 
between  x,  U,  and  W  by  using  the  expression  for  a z  in  Eq.  (1),  namely 


X  IdU  +  dv\ 
X+2G  \3x  dy  I 


(35) 


The  relationships  of  Eq.  (30)  between  stresses  and  strains  in  the  x  and  v  directions  are  also 
valid  here.  Since  the  parity  of  o  in  the  coordinate  z  is  even,  the  stress  integral  foxdz  ap¬ 
pears  in  the  equations,  unlike  the  case  of  bending  waves  where  the  moment  of  stress  is 
needed.  Thus  the  modulus  E/(l-v2)  in  Eq.  (30)  determines  the  wave  propagation,  the 
wave  equation  is  of  second  order,  and  the  phase  speed  for  extensional  waves  cp  is  inde¬ 
pendent  of  wavelength.  It  is  given  by 


-  . 

p  Pd-*'2) 

Pictorially  the  waves  display  flat  cross  sections  moving  to  and  fro  in  unison,  as 
sketched  in  Fig.  6. 
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Fig.  6  ~  Thin-plate  extensional  wave 


PLATE  STRESS  EQUATIONS  OF  MOTION 

The  development  of  approximate  theories  for  wave  propagation  in  plates  is  effected 
mathematically  by  integrating  the  equations  of  motion  (Eq.  (3)).  This  was  mentioned  be¬ 
fore  in  connection  with  the  discussion  of  thin-plate  theory.  In  this  section  the  derivation 
of  the  integrated  equations  is  given  explicitly,  to  serve  as  a  basis  for  the  understanding  of 
the  existing  thick-plate  theories  and  to  give  the  background  for  evaluation  of  correction 
factors  in  these  theories.  These  integrated  equations  are  often  called  plate  stress  equations. 
Such  an  integrated  set  of  equations  does  not  by  itself  constitute  a  new  theory  or  a  sim¬ 
plification.  Approximate  expressions  for  the  basic  variables  as  given  in  Eqs.  (27)  and 
possible  further  restrictions  are  discussed  in  connection  with  the  thin-plate  theories  are 
needed  to  make  the  mathematics  more  tractable. 

Like  any  function,  the  displacement  functions  u,  v,  w  may  be  written  as  a  sum  of 
an  odd  and  an  even  part  with  respect  to  the  coordinate  z.  Antisymmetric  waves  have  odd 
u,  v  and  even  w ;  symmetric  waves  have  even  u,  v  and  odd  w.  The  parity  of  the  deriva¬ 
tives  of  the  stresses  in  Eq.  (3),  in  view  of  Eqs.  '(1)  and  (2),  is  such  that  upon  integra¬ 
tion  across  the  thickness  direction  of  the  plate  the  first  two  equations  contain  functions 
belonging  only  to  the  antisymmetric  type,  whereas  the  third  contains  functions  belonging 
only  to  the  symmetric  type.  To  obtain  a  complete  set  of  equations  for  either  case  re¬ 
quires  forming  integrated  equations  from  the  moment  relationships,  which  follow  from 
Eq.  (3)  by  multiplying  every  equation  by  z.  Thus  two  independent  sets  of  equations  are 
derived  (one  for  each  type  of  wave),  as  is  shown  more  extensively  below.  If  one  assumes 
that  the  plate  is  not  subject  to  external  normal  and  shear  stresses,  the  two  types  of  waves 
can  occur  independently.  Nonhomogeneous  boundary  conditions  cause  a  mixture  of  the 
two  types  of  waves.  In  this  report  only  homogeneous  boundary  conditions  are  assumed. 

For  antisymmetric  waves  the  moments  are  readily  interpreted  as  torques  of  the 
stresses  and  angular  momenta  of  the  displacements.  Such  an  identification  is  not  obvious 
in  the  case  of  symmetric  waves,  and  the  moments  appear  only  as  mathematical  moments 
in  the  ^  variable  of  the  various  functions. 


Antisymmetric  Waves 

In  this  case  the  parity  of  the  functions  w,  oxz,  and  a  are  even  in  ^  while  u,  v, 
o x ,  ny ,  and  »  are  odd  in  ^.  Therefore,  the  pertinent  plate  stress  equations  are  those 

where  one  takes  the  moment  equation  following  from  the  first  two  parts  of  Eq.  (3)  and 
the  force  equation  following  from  the  third  part  of  Eq.  (3).  As  a  result,  one  finds  that 
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and 


3Af, 

d* 


dMxy  „ 
by  *x 


bMvx  bMv  „  b2Lv 

— i^-+ — 2-  -  Q  =  a - -2- 

bx  by  y  bt 2 


bQx  bQ 

_ + _ y 

dx  by 


(37) 


where  one  introduces  two  bending  moments  Mx ,  My ,  one  twisting  moment  A/xy ,  two 

transverse  shear  forces  Qx,  Qy,  one  integrated  vertical  displacement  Tz,  and  two  moments 

of  displacement  L  ,  L  .  These  symbols  are  defined  by  the  expressions 
*  y 


and 


Symmetric  Waves 

In  this  case  the  parity  of  the  functions  with  respect  to  z  is  u,  v,  ox,  oy,  oz,  axy  even 
and  w,  oxz,  oyz  odd.  Therefore  the  pertinent  plate  stress  equations  are  those  obtained  by 
the  direct  integration  of  the  first  two  parts  of  Eq.  (3)  and  from  the  moment  equation 
connected  with  the  third  part  of  Eq.  (3).  One  finds  that 


NRL  REPORT  8372 


and 


dNx 

dNxy 

d2T 

=  n  x 

dx 

dy 

P  dt* 

dN 

d2T 

dx 

dy 

p  dt 2 

(39) 


in  terms  of  three  normal  stress  integrals  Nx,  Ny,  Nz,  one  shear  stress  integral  Nxy<  and 
two  shear  stress  moments  Rx,  Ry.  Further,  one  needs  two  integrated  displacements 
Tx,  Ty  and  one  displacement  moment  Lz : 


Ux  = 


N*>m 


R _  = 


r  = 


i  to l 


and 


+d 

zwdz. 


THICK-PLATE  THEORY 


The  set  of  integrated  Eqs.  (37)  and  (39)  for  wave  propagation  in  a  plate,  with  the 
definitions  of  Eqs.  (38)  and  (40),  do  not  constitute  by  themselves  a  simplification  of  the 
problem.  Only  if  one  makes  assumptions  concerning  the  displacement  components  u,  v,  w 
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that  are  plausible  for  a  plate  can  one  reduce  these  integrated  equations  to  a  form  that  is 
readily  solved.  Such  simplifying  assumptions  were  already  given  in  Eq.  (27).  Thin-plate 
theory  is  characterized  by  further  restrictions  relating  the  functions  4>x,  <t>y,  X  to  U,  V, 

W  in  Eqs.  (29)  and  (35).  Models  of  wave  propagation  in  plates  based  on  £he  linear  ex¬ 
pansion  of  u,  v,  w  in  terms  of  z,  Eq.  (27),  without  further  restrictions,  are  indicated  by 
thick-plate  theory.  Of  course  it  is  quite  conceivable  that  one  would  try  to  retain  terms  of 
higher  order  than  the  first  in  the  series  expansion  of  u,  v,  w  to  improve  the  accuracy  of 
the  approximation.  Such  higher  order  models  have  so  far  not  come  to  the  attention  of 
the  author  (see  also  Appendix  C). 

Timoshenko-Mindlin  Plate  Theory 

The  dispersion  relation  for  antisymmetric  waves  in  thin  plates,  Eq.  (33),  has  an  un¬ 
realistic  property  in  the  sense  that  the  phase  speed  increases  beyond  bounds  as  fed  becomes 
very  large,  i.e.,  when  the  wavelength  becomes  small  compared  with  the  plate  thickness. 

The  reason  for  this  is  that  the  entire  plate  cross  section  in  the  thin-plate  model  is  sup¬ 
posed  to  partake  in  a  bending  motion.  As  a  consequence,  the  bending  stiffness  D  in 
Eq.  (31)  is  proportional  to  the  third  power  of  the  plate  thickness.  The  plate  becomes  so 
rigid  when  d  increases  that  a  disturbance  is  propagated  instantaneously,  which  makes  this 
model  not  acceptable  on  physical  grounds. 

Timoshenko  found  a  means  of  circumventing  a  similar  nonphysical  behavior  of  the 
phas-1  speed  in  the  case  of  bars.  The  same  basic  idea  of  Timoshenko  was  applied  to  plates 
by  Mindlin  [1].  Physically,  this  idea  amounts  to  restoring  rotatory  inertia  to  its  role  in 
the  wave  propagation  in  the  plate  and  to  dropping  the  requirement  that  cross  sections 
originally  perpendicular  to  the  neutral  plane  remain  so  during  bending.  Both  effects, 
rotatory  inertia  and  shear  strain,  were  ignored  in  the  derivation  of  bending  waves  in  a 
thin  plate.  Mathematically  the  proper  equations  are  obtained  by  inserting  the  set  of 
Eq.  (28)  into  the  force,  moment,  and  displacement  integrals  without  the  restrictions  of 
Eq.  (29).  This,  with  Eq.  (30),  results  in 


and 

Tz  =  2 dW. 
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The  symbol  (>'  plays  the  role  of  an  “effective  shear  modulus”  and  accounts  for  the  fact 
that  the  linear  expansion  in  2  of  the  u,  w  functions  is  only  an  approximation.  Without 
its  introduction  the  dispersion  relation,  although  bounded,  would  not  approach  a  value 
plausible  in  the  light  of  exact  elasticity  theory,  namely  the  Rayleigh  wave  speed  cJ(  If 
one  writes 

O'  =  k]G,  (42) 

the  effective  shear  modulus  is  represented  by  a. dimensionless  constant  k,  . 


Further  development  of  the  Timoshenko-Mindlin  theory  consists  in  making  the  in¬ 
sertion  of  F.q.  (41)  into  F,q.  (37),  differentiating  the  first  equation  of  Kqs.  (37)  with  respect 
to  .v  and  the  second  with  respect  to  y,  and  then  adding  the  results  (see  Ref.  8  for  details). 

<1  (p  (10 

There  results  two  equations  in  the  quantities  <l>  =  - —  +  — *-  and  W  . 

<)x  3v 


D  2<I>  -  2k)Gc1<1>  -  2k |  Gd  2  W  = 


2 pd*  i) 2  <l‘ 
3  bC- 


and 


(43) 


2n]Gd\  2W  +  <!>] 


2 pd 


d2W 

bt2 


By  assuming  a  wave  solution  for  T  and  IF  of  the  form  exp[/(cut-fe.v  I] ,  one  t intis  tin 
dispersion  relationship 


~(ltd)2c2  -  —(ltd)2  c2  -  k2c2 
,3  3  p 


which  can  be  written  as 

-  *  (kd)2( y2  -  7p)(72  -  k]  )  =  k2\V  ■ 


=  0.  (44) 


(45) 


The  symbol  y  is  used  for  the  ratio  c/cs .  The  subscripts  for  7  follow  the  subscripts  for 
the  phase  speed  c. 

This  equation  is  quadratic  in  the  variable  y2 .  The  smaller  root  can  be  identified 
with  the  zero-order  Lamb  wave.  At  low  frequencies  (i.e.,  ltd  close  to  zero)  this  root  gives 
a  dispersion  relationship 


-(fed)2  7^,  or 


c  =  ch 


(46) 
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which  is  the  dispersion  relation  for  bending  waves  in  a  thin  plate  (Eq.  (33)),  as  one  would 
expect.  This  same  root  has  a  horizontal  asymptote  given  by 

y2  =  .  (47) 

This  result  offers  Mindlin  the  means  to  fix  the  value  of  k 2 ,  since  it  is  physically  plausible 
that  the  thick-plate  phase  velocity  would  approach  the  phase  velocity  of  Rayleigh  waves, 
which  are  waves  occurring  at  the  surface  of  a  semi-infinite  solid.  This  choice  of  does 
not  reproduce  the  fact  that  this  correction  factor  should  be  frequency  dependent.  As  a 
consequence,  it  might  be  expected  that  the  dispersion  relation  with  this  value  of  n]  will 
not  be  very  close  to  the  exact  relation  (Eq.  (23))  at  intermediate  frequencies;  but,  indeed, 
the  agreement  is  quite  good  (see  Tables  1—3). 


Table  1  —  Dispersion  Relation  for  Antisymmetric 
Waves  According  to  Various  Theories 


kd 

y,  I* 

y,  //t 

7.  HI* 

0.3 

0.2319 

0.2384 

0.2387 

0.9 

0.5544 

0.5449 

0.5478 

2.3 

0.7963 

0.7696 

0.7766 

3.9 

0.8607 

0.8294 

0.8388 

5.9 

0.8795 

0.8517 

0.8626 

7.5 

0.8827 

0.8602 

0.8703 

8.7 

0.8834 

0.8648 

0.8737 

10.1 

0.8836 

0.8690 

0.8762 

*/  —  Exact  theory. 

^7/  —  Thick-plate  correction  factor  proposed  in  this 
report. 

Ill  —  Thick-plate  correction  factor  proposed  by 
Mindlin  (1  ] . 

Note:  v  =  0,05 

yR  =  0.8837 
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Table  2  —  Dispersion  Relation  for  Antisymmetric 
Waves  According  to  Various  Theories 


kd 

7,  /* 

7,  //* 

7.  ///' 

0.3 

0.2721 

0.2751 

0.2763 

0.7 

0.5338 

0.5245 

0.5322 

1.1 

0.6801 

0.6614 

0.6758 

1.5 

0.7643 

0.7384 

0.7575 

2.5 

0.8602 

0.8256 

0.8496 

4.5 

0.9122 

0.8773 

0.9003 

6.5 

0.9240 

0.8948 

0.9142 

8.5 

0.9269 

0.9041 

0.9197 

10.1 

0.9276 

0.9091 

0.9220 

*1  —  Exact  theory. 

^ II  —  Thick-plate  correction  factor  proposed  in  this 
report 

111  —  Thick  plate  correction  factor  proposed  by 
Mindlin  [  1  ] . 

Note:  V  =  0.3028 
yR  =  0.9278 
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Table  3  —  Dispersion  Relation  for  Antisymmetric 
Waves  According  to  Various  Theories 


fed 

7,  I* 

7.  II* 

7,  HI* 

0.3 

0.4827 

0.3196 

0.3222 

0.5 

0.5998 

0.4759 

0.4841 

0.9 

0.7424 

0.6642 

0.6854 

1.3 

0.8180 

0.7589 

0.7893 

1.9 

0.8771 

0.8279 

0.8654 

2.5 

0.9064 

0.8619 

0.9016 

3.3 

0.9269 

0.8866 

0.9256 

4.3 

0.9397 

0.9040 

0.9400 

5.3 

0.9465 

0.9147 

0.9474 

6.9 

0.9517 

0.9254 

0.9533 

8.9 

0.9541 

0.9332 

0.9570 

9.9 

0.9546 

0.9360 

0.9580 

*/  —  Exact  theory. 

T/J  —  Thick-Plate  correction  factor  proposed  by  this 
report. 

Ill  —  Thick-Plate  correction  factor  proposed  by 
Mindlin  [X], 

Note:  v  =0.5 

yR  =  0.9553 

Mindlin  [1]  does  not  consider  the  larger  root  of  Eq.  (45).  In  the  section  on  gen¬ 
eralized  correction  factors  it  is  shown  that  this  root  can  be  identified  with  the  first-order 
mode  in  antisymmetric  Lamb  waves. 


Symmetric  Waves  in  Thick  Plates 

Although  the  symmetric  waves  in  a  thin  plate  do  not  have  the  singularity  for  fed  =  0 
of  the  bending  waves,  the  constant  phase  speed  cp,  given  by 
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E  'A 
p(l-p2) 

does  not  lead  to  the  expected  asymptotic  behavior,  for  high  frequencies,  namely  the 
Rayleigh  wave  speed.  This  suggests  that  a  correction  should  be  developed  for  symmetric 
waves  in  a  thin  plate  parallel  to  the  correction  introduced  by  Mindlin  for  antisymmetric 
waves  in  a  thin  plate. 

Such  an  effort  was  indeed  carried  out  by  Kane  and  Mindlin  [2]  but  from  a  dif¬ 
ferent  standpoint.  The  authors  point  out  in  their  paper  that,  on  the  basis  of  thin-plate 
theory,  the  characteristic  vibrations  of  a  finite  plate  are  determined  by  the  lateral  dimen¬ 
sions  of  the  plate.  If  the  diameter  of  a  circular  plate  becomes  small  compared  with  the 
thickness,  the  thin-plate  theory  does  not  give  the  proper  low-frequency  modes,  and  one 
has  to  introduce  thickness  vibrations  in  order  to  obtain  realistic  values  for  the  lowest 
modes  of  vibration.  This  is  a  different  line  of  inquiry  than  the  one  followed  in  this  study, 
where  the  mechanism  of  a  thin-plate  wave  is  considered  in  the  transition  from  a  true 
volume  phenomenon  toward  the  surface  waves  of  Rayleigh.  Therefore  a  further  discus¬ 
sion  of  the  Kane-Mindlin  method  will  not  be  offered  here.  Instead,  a  development  analo¬ 
gous  to  the  one  above  for  the  Timoshenko-Mindlin  theory  will  be  presented,  leading  to 
the  formulation  of  another  effective  shear  modulus  —  in  this  case  for  symmetric  waves 
in  a  thick  plate. 

One  is  guided  to  this  formulation  by  the  parallel  case  of  antisymmetric  waves,  as 
developed  in  the  section  on  Timoshenko-Mindlin  plate  theory.  One  enters  the  series 
expansion  (Eq.  (34))  for  the  displacement  components  in  terms  of  the  coordinate  z  into 
the  integrals  in  Eqs.  (40)  that  represent  integrals  of  stresses,  moments  of  stresses,  and 
components  of  displacement.  Of  course,  the  thin-plate  approximation  (Eq.  (35))  is  not 
applied  here.  The  stress  components  are  expressed  in  terms  of  derivatives  of  displace¬ 
ment  components  by  means  of  Eqs.  (1).  The  result  is  a  set  of  equations  parallel  to 
Eqs.  (41)  for  antisymmetric  waves: 

[{2  ♦(JA_)(|K.x) 

a-  m 

R  =  —  d3G"-^&--  R  =  —  d3G"~* 
x  3  bx  ’  y  3  by 

Lz  =  -|d3x;  Tx=  2d 17;  Ty  =  2dV. 
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In  searching  for  a  correction  to  extensional  waves  in  thin  plates,  one  would  be  inclined 
to  follow  the  successful  example  of  the  introduction  of  transverse  shear  and  rotatory 
inertia  by  Mindlin  for  antisymmetric  waves,  although  admittedly  the  latter  factor  had 
much  less  influence  on  the  dispersion  relation  than  the  former.  The  major  improvement 
proposed  by  Mindlin  is  the  introduction  of  an  effective  shear  modulus  G'  in  the  shear 
stress  integrals  Qx  and  Qy.  The  corresponding  integrals  in  symmetric  waves  are  integrals 
of  moments  of  shear  stress,  Rx  and  R  .  It  appears  plausible  to  assume  that  a  similar 
improvement  would  result  from  introducing  another  effective  shear  modulus  G"  for 
symmetric  waves  in  the  expression  of  these  integrals,  as  is  shown  in  Eq.  (49).  The  fol¬ 
lowing  analysis  shows  the  validity  of  this  choice.  A  nondimensional  correction  factor  is 
introduced  to  represent  this  effective  modulus  by 

G"  =  k\G.  (50) 

Equations  (49)  and  (50)  are  inserted  into  the  equations  of  motion  (Eqs.  (39)).  If  one 
differentiates  the  first  equation  of  motion  with  respect  to  x,  and  the  second  with  respect 
to  y,  and  adds  the  resulting  equations,  one  obtains 


(X+2G)v2'!'  +  Xv2x  =  P 


a2* 

bt2 


and  the  third  equation  of  motion  becomes 

-XT'  +~-«2Gd2v2x  -  (X+2G)x  =  > 


(51) 


(52) 


where 

dx  dy 

Assuming  that  T'  and  x  are  represented  by  a  straight-crested  wave  in  the  form 
exp[i(cot-fex)] ,  one  obtains  the  dispersion  relation 

I  X+2G-pc2  X  1 


X 


22G(kd)2  --i(pc2)(fed)2  +  X  +  2G 
o  o 


(53) 


If  one  introduces  the  various  wave  speeds  corresponding  to  the  elastic  moduli  in  this 
equation  and  represents  the  wave  speeds  in  dimensionless  form  by  y  =  c/cg,  then  the 
dispersion  relation  appears  in  the  alternative  form, 


^  -  y2  K  - 2 

y\  -  2  ^K\{kd)2  -\y2{kd)7  +  y2 

6  6 


(54) 
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This  dispersion  relation  is  quadratic  in  the  variable  y2 ,  as  in  the  case  of  antisymmetric 
waves.  The  smaller  root  can  be  identified  with  the  zero-order  symmetric  Lamb  wave. 
At  low  frequencies  (i.e.,  kd  close  to  zero),  this  root  gives  for  the  dimensionless  wave- 
number 


y2  *  4(1  -  1/7*)  =  7*.  (55) 

This  is  the  phase  speed  for  extensional  waves,  as  derived  in  the  theory  of  symmetric 
waves  in  thin  plates.  This  same  root  has  a  horizontal  asymptote,  given  by 

72=k|.  (56) 

Thus  by  a  proper  choice  or  definition  of  *c2  (as  discussed  below),  one  can  ensure 
that  the  phase  speed  in  thick-plate  theory  (for  symmetric  waves)  also  approaches  the 
phase  speed  of  Rayleigh  waves.  One  sees,  therefore,  that  the  factor  plays  a  similar 
role  in  fitting  the  high-frequency  end  of  the  dispersion  relation  for  symmetric  waves,  as 
the  factor  Kt  did  in  the  case  of  antisymmetric  waves.  Both  are  correction  factors  to  the 
shear  modulus  G. 

Just  as  in  the  case  of  antisymmetric  waves,  one  can  identify  the  second  root  of  the 
dispersion  relation  (Eq.  (54))  with  the  first-order  mode  of  symmetric  Lamb  waves.  This 
is  discussed  more  extensively  below. 


EFFECTIVE  SHEAR  MODULUS  IN  ANTISYMMETRIC  WAVES 

The  coefficient  k?  defining  the  effective  shear  modulus  with  respect  to  the  actual 
shear  modulus  in  antisymmetric  waves  (Eq.  (42))  was  fixed  by  Mindlin  [1]  at  a  constant 
value  such  that  the  phase  speed  asymptotically  approaches  the  Rayleigh  wave  speed. 

It  is  possible  to  derive  the  value  of  k}  directly  from  theory  by  comparing  the  ex¬ 
pression  for  the  transverse  shear  force  Qx  from  the  plate  stress  equations  (Eqs.  (37)  and 
(38))  by  utilizing  Eq.  (1), 


Q 


X 


(57) 


with  the  corresponding  expression  in  the  thick-plate  equations  (Eq.  (41))  and  using 
Eq.  (42), 


Qx  = 


2k}  Gd 


(dW 
\  dx 


) 


(58) 


This  shows  that  the  correction  factor  k2  should  be  defined  by 
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(59) 


The  numerator  in  this  definition  can  be  obtained  from  exact  elasticity  theory.  The  de¬ 
nominator  contains  the  quantities  that  are  the  field  variables  in  thick-plate  theories.  To 
relate  them  to  the  exact  theory,  consider  the  defining  equations  (Eq.  (28)) 

u(x,2,t)  =  z<Px(x,t) 


and  (60) 

w(x,z,t)  =  W'fx.t). 

It  follows  that  the  derivative  of  the  vertical  average  displacement  is  given  as 


(61) 


The  angle  <Px  is  the  average  angle  of  rotation  of  a  plate  cross  section.  It  appears  plausible 
to  relate  this  average  angle  to  the  displacement  component  u  from  elasticity  theory 
through  the  expression 


This  choice  for  <S>X  leads  to  the  desired  asymptotic  behavior  of  the  phase  speed. 


(62) 


The  detailed  formulae  for  the  quantities  in  the  Eqs.  (59),  (61),  and  (62),  in  terms 
of  the  parameters  of  the  Lamb  waves,  are  given  in  Appendix  A.  From  Eqs.  (A8),  (A9), 
and  (A10)  one  can  infer  the  asymptotic  values  for  the  parts  of  Eq.  (59),  considering 
that  for  large  argument  the  hyperbolic  functions  approach  the  exponential  function. 

Since  for  the  hyperbolic  integral,  shi(x),  defined  as  shi(x)  =  /*  sinh(f)/fdt,  ap¬ 
proaches  the  limit  exp  (x/x)  (see  Ref.  9),  the  value  of  ax  in  the  denominator  of  Eq.  (59) 
can  be  ignored  in  comparison  with  the  value  of  dW/dx,  in  the  high-frequency  limit. 


Then  for  one  has 


i 
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Ca\  (fe2+s2(s-q) 
d *}  qs 


exp(sd) 


(63) 


and 


(C\  (fe2s-2fe2q+s3 ) 
\c/2/  2  qs 


exp(sd). 


(64) 


As  a  consequence,  the  asymptotic  value  of  the  correction  factor  K2  is  given  by 


K. 


1 

1 


2 (s-q)(k2+s2 ) 
fe2s  +  s3  -  2 k2q 


for  fed 


OO 


(65) 


On  the  other  hand,  the  dispersion  relation  for  Rayleigh  waves  derived  in  exact 
elasticity  theory,  Eq.  (25),  can  be  transformed  into  the  form 

2(s-<?)(fe2+s2)  fe2-s2 

fe2s  +  s3  -  2k2 q  ~  k 2  (  b) 

by  algebraic  manipulation.  Since  (fe2-s2)/fe2  =  y2  and  the  value  of  y2  for  Rayleigh  waves 
is  the  relative  phase  speed  yR,  one  see'  by  comparing  Eqs.  (65)  and  (66)  that  in  the  limit 
of  high  frequency, 


-*■  Tfi  ,  for  (fed)  -*  oo.  (67) 

This  is  exactly  the  proper  high-frequency  behavior  of  the  correction  factor.  Here,  how¬ 
ever,  the  required  high-frequency  behavior  is  not  imposed  as  in  the  Timoshenko-Mindlin 
theory,  but  follows  naturally  from  the  definition  of  k2  as  given  by  Eq.  (59),  in  terms 
of  results  from  exact  elasticity  theory. 

The  results  of  calculating  the  correction  factor  k,  according  to  Eqs.  (59),  (61),  and 
(62)  are  presented  in  Fig.  7.  Three  values  of  Poisson’s  ratio  were  chosen  to  show  the 
variation  of  k  ,  as  a  function  of  this  parameter.  The  asymptotic  value,  equal  to  the 
relative  Rayleigh  wave  speed  yR ,  is  indicated  by  an  arrow. 

One  can  calculate  the  relative  wave  speed  in  the  thick-plate  approximation  as  a 
function  of  the  relative  wavenumber  fed  from  Eq.  (45).  If  one  chooses  the  constant  value 
*i  =  yR  proposed  by  Mindlin,  the  result  is  quite  close  to  the  exact  dispersion  relation, 

Eq.  (23).  If  one  chooses  for  the  correction  factor  k,  in  the  dispersion  relation,  Eq.  (45), 
the  frequency-dependent  values  from  Fig.  7,  the  agreement  with  exact  theory  is  less 
favorable,  except  in  a  few  cases.  This  comparison  is  shown  in  Tables  1,  2,  and  3. 

i 
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Fig.  7  —  Correction  factor  K  i  for  effective  shear  modulus  in  antisymmetric  waves;  (1)  e  =  0.05, 
(2)  v  =  0.3028,  (3)  v  =  0.5.  The  arrow  -*  indicates  the  asymptote. 


EFFECTIVE  SHEAR  MODULUS  IN  SYMMETRIC  WAVES 

It  is  possible  to  derive  an  expression  for  the  effective  shear  modulus  in  symmetric 
waves  in  a  way  similar  to  the  case  of  antisymmetric  waves.  Compare  the  expression  for 
the  stress  moment  integral  Rx  as  given  in  Eq.  (40),  utilizing  Eq.  (1), 


R 


X 


(68) 


with  its  expression  as  given  in  thick-plate  theory,  Eq.  (49),  utilizing  Eq.  (50), 

it  ,69) 

This  shows  that  the  correction  factor  nj  should  be  defined  by 
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«22 


(70) 


The  numerator  in  this  definition  can  be  obtained  directly  from  exact  elasticity 
theory.  The  denominator  contains  the  variable  x  from  thick-plate  theory.  To  relate  this 
quantity  to  the  exact  theory,  consider  the  defining  Eqs.  (34), 


u(*,^,0  =  U{x,t) 


and 


(71) 


Ul(x,Z,t)  =  2\(X,t). 

There  are  several  ways  to  define  the  displacement  coefficient  X-  The  choice 


J>L-  (72, 

dx  UJ  \dxj 
0 

appears  plausible  and,  moreover,  leads  to  the  desired  asymptotic  behavior  of  the  wave 
speed. 

The  detailed  formulae  for  the  quantities  in  Eqs.  (70)  and  (72),  in  terms  of  the 
parameters  of  the  Lamb  waves,  are  given  in  Appendix  A. 

The  numerator  of  Eq.  (70)  is  given  by  Eq.  ( A 1 9 ) ;  and  the  denominator,  as  further 
defined  by  Eq.  (72),  is  given  by  Eq.  (A20).  Equations  (A19)  and  (A20)  are  obtained  by 
integration  by  parts  and,  as  a  consequence,  contain  two  parts  such  that  the  second  part  is 
smaller  by  one  order  in  kd  than  the  first  part.  Therefore,  the  second  part  is  negligible 
in  comparison  with  the  first  part  in  the  limit  of  large  hd.  Since,  moreover,  the  hyperbolic 
functions  approach  the  exponential  function  for  large  kd,  one  sees  that  these  first  parts 
become  equal  to  the  formulae  (Eqs.  (A8)  and  (A9))  for  antisymmetric  waves,  respectively. 
Thus,  the  proof  given  before  for  the  asymptotic  value  of  the  correction  factor  in  antisym¬ 
metric  waves  is  the  same  for  the  correction  factor  in  symmetric  waves.  This  means  that  also 
in  the  latter  case  kJ  approaches  the  desired  value  7 ^  for  kd  -*  However,  the  behavior 
of  K 2  for  intermediate  and  lower  kd  is  unsatisfactory,  in  the  sense  that  it  does  not  stay 
positive.  This  difference  in  behavior  as  compared  with  the  case  of  antisymmetric  waves 
is  due  to  the  factors  (ltd)2  +  (xd):  and  (s</):.  The  dispersion  relation  for  symmetric  waves 
is  such  that  these  factors  become  negative  for  low  values  of  kd.  1'he  conclusion,  there¬ 
fore,  has  to  be  that  the  definition  of  k{  according  to  Eq.  (70)  combined  with  Eq.  (72i 
does  provide  the  correct  high-frequency  value,  but  it  cannot  be  used  for  lower  frequencies. 
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Efforts  to  improve  this  situation  should  he  considered  against  the  background  of  general¬ 
ized  correction  factors.  This  is  discussed  in  Appendix  B. 


The  dispersion  relation  according  to  thick-plate  theory  for  symmetric  waves,  Eq.  (54), 
was  computed  using  a  constant  value  for  k  equal  to  yR.  The  results  are  given  in  Tables 
4,  5,  and  6  and  are  compared  with  the  dispersion  relation  for  symmetric  Lamb  waves, 

Eq.  (24).  One  can  observe  that  the  agreement  is  not  as  good  as  the  corresponding  case 
of  antisymmetric  waves,  Tables  1  through  3. 


Table  4  —  Dispersion  Relation  for 
Symmetric  Waves  According  to 
Various  Theories 


kd 

/* 

7/t 

0.1 

1.4509 

1.4509 

0.5 

1.4507 

1.4508 

0.9 

1.4501 

1.4505 

1.3 

1.4473 

1.4499 

1.7 

1.3413 

1.4479 

2.1 

1.1569 

1.4347 

2.5 

1.0504 

1.3329 

2.9 

0.9878 

1.2340 

3.3 

0.9500 

1.1653 

3.7 

0.9265 

1.1138 

4.5 

0.9023 

1.0449 

5.7 

0.8893 

0.9873 

6.9 

0.8855 

0.9557 

8.1 

0.8843 

0.9365 

9.7 

0.8838 

0.9208 

13.3 

— 

0.9036 

16.5 

— 

0.8967 

19.7 

— 

0.8928 

*/  —  Exact  theory. 

^ II  —  Thick-plate  theory,  a<  wording  to 
this  report,  with  constant  factor 

Note:  V  *  0.05 

yH  =  0.8837 
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Table  5  —  Dispersion  Relation  for 
Symmetric  Waves  According  to 
Various  Theories 


kd 

/* 

//t 

0.1 

1.6887 

1.6933 

0.5 

1.6626 

1.6841 

0.9 

1.5932 

1.6601 

1.3 

1.4460 

1.6158 

1.7 

1.2750 

1.5465 

2.1 

1.1507 

1.4586 

2.5 

1.0717 

1.3698 

2.9 

1.0222 

1.2928 

3.3 

0.9909 

1.2305 

3.7 

0.9707 

1.1810 

4.5 

0.9485 

1.1106 

5.7 

0.9353 

1.0480 

6.9 

0.9307 

1.0124 

8.1 

0.9290 

0.9904 

9.7 

0.9282 

0.9722 

13.3 

— 

0.9519 

16.5 

— 

0.9437 

19.7 

— 

0.9391 

*/  —  Exact  theory. 

^7/  —  Thick-plate  theory,  according  to 

this  report,  with  constant  correction 
factor  k2  ~  yR- 

Note:  v  =  0.3028 
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Table  6  —  Dispersion  Relation  for 
Symmetric  Waves  According  to 
Various  Theories 


fed 

I* 

//+ 

0.0 

— 

2.0000 

0.4 

1.9467 

1.9607 

0.8 

1.8017 

1.8600 

1.2 

1.5986 

1.7329 

1.6 

1.4012 

1.6077 

2.0 

1.2519 

1.4978 

2.4 

1.1519 

1.4067 

2.8 

1.0872 

1.3330 

3.2 

1.0455 

1.2738 

3.6 

1.0181 

1.2264 

4.4 

0.9874 

1.1570 

5.6 

0.9683 

1.0930 

6.8 

0.9610 

1.0552 

8.0 

0.9579 

1.0314 

9.6 

0.9563 

1.0115 

13.2 

— 

0.9799 

16.4 

— 

0.9747 

19.6 

JZ _ 

0.9747 

*1  —  Exact  theory. 

^ II  —  Thick-plate  theory,  according  to  this 
report,  with  constant  factor  K2  =  ~tR  ■ 

Note:  V  =0.5 

yR  =  0.9553 

CONCLUSIONS  AND  PLANS  FOR  FURTHER  WORK 

The  main  result  of  the  study  herein  reported  is  the  determination  of  a  correction 
factor  for  effective  shear  modulus  in  antisymmetric  and  symmetric  waves  in  thick  plates 
by  comparison  with  elasticity  theory.  This  may  be  contrasted  with  Mindlin’s  (1]  ap¬ 
proach  for  antisymmetric  waves,  whereby  the  correction  factor  is  fit  to  represent  the 
high-frequency  limit  of  the  dispersion  relation  in  such  a  way  that  the  phase  speed  is 
equal  to  the  Rayleigh  wave  speed.  The  correction  factor  proposed  in  this  report  ac¬ 
complishes  this  without  recourse  to  arguments  outside  the  theory  proper,  as  in  Mindlin  [1] 
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it  can  be  shown  analytically  that  the  phase  speed  based  on  this  correction  factor  approaches 
the  Rayleigh  wave  speed  asymptotically.  This  obtains  both  for  antisymmetric  and  sym¬ 
metric  waves. 

This  approach  can  be  extended  to  the  correction  of  approximate  terms  in  the  thick- 
plate  theories  other  than  the  term  represei.ting  transversal  shear.  Certain  difficulties  with 
singularities  arise  that  have  not  been  satisfactorily  resolved  thus  far.  The  method  promises 
to  more  accurately  identify  the  larger  root  in  the  quadratic  dispersion  relation  from  thick- 
plate  theory  with  the  first  higher  mode  of  antisymmetric  and  symmetric  Lamb  waves 

In  addition  to  further  study  of  this  generalization  of  the  methods  in  this  report,  the 
next  step  contemplated  at  this  time  is  the  evaluation  of  correction  factors  for  effective 
shear  modulus  in  cases  where  the  plate  is  subject  to  other  than  homogeneous  boundary 
conditions:  a  plate  loaded  on  both  sides  or  one  side  by  a  fluid,  a  composite  plate,  and  a 
composite  plate  under  fluid  loading. 
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FORMULAE  FOR  FIELD  VARIABLES  FROM  ELASTICITY  THEORY 


The  symbols  used  in  this  appendix  are  defined  in  the  text  of  the  report  and  also  can 
be  found  alphabetically  arranged  in  the  list  of  symbols  on  pages  v-viii. 

The  dependence  on  time  and  at -coordinate  of  all  variables  occurs  through  a  factor 
exp[i(cof-/iA:)]  that  is  not  repeated  in  the  following  formulae.  The  amplitudes  CQ  and 
Ds  have  the  dimension  of  length  squared; 

Q 1  ~  k2  [1  -(c/cd)2  ]  s2 


s  2p(l+f) 

Antisymmetric  Waves 

For  antisymmetric  waves  one  uses  the  part  of  the  potential  0  that  has  odd  parity  in 
z,  and  the  part  of  the  potential  \p  that  has  even  parity  in  z : 

0  =  jBq  sinh(qz)  0  =  CQ  cosh(sz)  .  (A2) 

The  value  of  c  for  given  kd  follows  from  the  dispersion  relation  for  homogeneous 
boundary  conditions  (Eq.  (23)); 

(k2  +s2  )2sinh(qd)cosh(scf)-4/?2<?scosh(<jd)sinh(sd)  =  0.  (A3) 

From  Eq.  (22),  use  the  relation 

Ra2ikq  cosh(q d)+Ca(k2  +sJ  )cosh(sd)  =  0  (A4) 


=  k2  [  l-(c/cs)2  ] 

= _ £(!-»>) 

p(l  ^)(l-2e) 


(Al) 


to  derive 

u_l/30_30\ 
d  d  \dx  bz! 

(A5) 

j\(kd)2 +(sd)2  lcosh(sd)sinh(qz)-2(qd )(sd)cosh(qd)sinh(sz)  \ 
d2  ^  2( qd  )cosh(qd )  I 
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bw  32  0  +  32  0 
bx  bxbz  3x2 


Ca  ([(fed)2 +(sd)2  ]cosh(sd)cosh(q2)-2(fed)2cosh(qd)cosh(s2) 
d 2  1  2cosh(gd) 


e 


zx 


[(fcd)2+(sd)2] 


cosh(sd)cosh(<j2)-cosh(<jd)cosh(s2) 

cosh(gd) 


(A6) 


( A7) 


1/ 


.d2 


Is.  [(fed)2  +(sd)2 


(sd  )cosh  (sd  )sinh  (g  d  )-(gd  )cosh(  qd  )si  nh  (sd )  1 


(A8) 


(qd)(sd)cosh(qd) 


Cfl  ((sd)[(fed)2  +(sd)2  ]  cosh(sd)sinh(gd)-2(fed)2  (gd)sinh(sd)cosh(qd) 
d2  '  2(<jd)(sd)cosh(gd) 


(A9) 


and 


d2 

Ca  j[(kd)2  +(sd)2  ]  cosh(sd)shi(qd)-2(qd)(sd)cosh(gd)shi(sd) 
d2  '  2(qd)cosh(gd) 

The  hyperbolic  integral  shi(x)  in  Eq.  (A10)  is  defined  by 

shi(x)  =f  2f-dt 

and  is  computed  according  to  the  series  expansion 


(A10) 


(All) 


34 


NRL  REPORT  8372 


Symmetric  Waves 


,2n+l 


shi(x)  = 


'(2n+l)(2n+l)! 


(A12) 


For  symmetric  waves  one  uses  that  part  of  the  potential  <p  that  has  even  parity  in  z, 
and  the  part  of  the  potential  \p  that  has  odd  parity  in  z: 


<t>  -  A4cosh(gz);  0  =  Dssinh(sz) 

The  value  of  c  for  a  given  fed  follows  from  the  dispersion  relation  for  homogeneous 
boundary  conditions  (Eq.  (24)); 


(fc2+s2)2cosh(gd)sinh(sd)-4fe2gssinh(qd)cosh(sd)  =  0. 


From  Eq.  (22),  use  the  relation 

As2ikqsinh(qd)+D$(k'1  +s2  )sinh(sd)  =  0 

to  derive 


u  _  1/90  djp\ 

d  d  \9x  9z  / 

Dg  j  [(fed)2  +(sd)2 1  sinh(sd)cosh(qz)-2(<jd)(sd)sinh(gd)cosh(sz) 
d2  '  2(qd)sinh(qd) 

9 w  _  920  +  9 20 

9x  9x9  z  9x2 


Ds  f[(fed)2+(sd)2  ]  sinh(sd)sinh(qz)-2(fed)2sinh(qd)sinh(sz)( 
d2  '  2sinh(qd)  f 

sinh(sd)sinh(qz)-sinh(qd)sinh(sz) 
sinh(qd) 


[(fed)2+(sd)2] 


(A14) 


(A15) 


(A16) 


(A17) 


(A18) 


i  r  d 

z€zxdz  =  -jt  f^)2+  (sdf 


(A19) 


X 


(sd)cosh(qd)sinh(sd)-(qd)sinh(qd)cosh(sd) 
(qd)(sd  )sinh(qd) 


(qd)2-(sd)2 

(qd)2(sd)2 


sinh(sd) 
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\f  7'— W  - 

7  \dx  I 


rA 


[(W)-j  (s f# ) -  |  (sd)<:osh(<7d)sinM.sd,-2(i?d)(/?d)2sinh(<2d)cosh(sd) 
2(<jd)(sd)sinh(<}d) 


uj  / 2( qd ) 2  (fed ) 2  -(sd ) 2  [  ( fed )  ~  +fsd)2 

'  2(/d)-'(sd)-’ 


•d):i| 


and 


I j  |  (fed)2  +(sd)2  ] sinh(a~d)shi(qd )-2(.sd i  sinh(<yd)shi(sd) 
■s  2sinh(r/d) 


(A20) 


(A21 ) 
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GENERALIZED  CORRECTION  FACTORS 


The  principle  by  which  the  correction  factors  in  the  foregoing  sections  were  derived 
can  be  generalized  to  cover  also  the  other  terms  in  the  equations  of  motion  for  waves 
in  thick-plate  theory,  Eqs.  (43)  (51),  and  (52).  Just  as  in  defining  the  correction  factor 
for  shear  modulus,  the  procedure  consists  in  comparing  the  integrals  defined  in  Eqs.  (38) 
and  (40)  with  the  thick-plate  approximations,  Eqs.  (41)  and  (49),  respectively.  This  shows 
that  the  approximations  Mx  and  Lx  have  essentially  the  same  correction  factor,  while 
the  correction  factor  for  T  is  equal  to  one.  Thus  only  one  additional  correction  factor  is 
needed  in  the  case  of  antisymmetric  waves  and  is  defined  by 


*1 


(Bl) 


This  correction  factor  is  the  same  for  the  term  representing  bending  stiffness  as  for  the 
rotary  inertia  in  Eq.  (43).  Applying  the  same  procedure  to  the  case  of  symmetric  waves 
leads  to  two  additional  correction  factors,  one  for  the  quantity  N  , 


and  one  for  the  quantity  N 


(B2) 


(B3) 


The  correction  factors  for  Tx  and  L__  are  equal  to  one,  provided  x  is  defined  as 
X  =  3/d3  £jd  zwdz.  The  factors  k*  and  k *  are  correction  factors  for  the  stiffness  modulus 
in  the  x  and  2  directions,  respectively.  The  subscripts  of  the  correction  factors  have  been 
chosen  such  that  those  belonging  to  antisymmetric  waves  are  odd  (K  ,  k  )  and  those  of 
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symmetric  waves  are  even  (#c2 »  *4.  K6)-  Inserting  the  correction  factors  k2,  k]  into  the 
equations  of  motion  for  antisymmetric  waves  (Eq.  (43))  leads  to  the  following  dispersion 
relation: 


i|ci(fed)V  ~^n](kd)2y2 


K 


2 

1 


75-«i 


=  0. 


(B4) 


The  smaller  root  of  this  quadratic  equation  is  identified  with  the  zero-order  antisym¬ 
metric  Lamb  wave.  The  larger  root  corresponds  to  the  first-order  Lamb  wave.  This  wave 
has  a  vertical  asymptote  if  one  represents  the  dispersion  relation  with  7  as  a  function  of 
kfd,  where  feg  is  the  wavenumber  of  the  shear  wave  corresponding  to  the  given  frequency 
By  use  of  the  relation  fegd  =  ykd,  Eq.  (B4)  is  transformed  into 


i»cS(fetd)V  -*i72  -r«3(fe4d)s72 


Ki72 


72-^Yi 


=  0. 


(B5) 


It  is  instructive  to  consider  the  behavior  of  this  dispersion  relation,  Eq.  (B5),  at  low  and 
at  high  frequencies  for  both  roots.  This  is  shown  in  Table  Bl.  For  comparison,  the  limits 
at  low  and  at  high  frequencies  of  the  dispersion  relation  from  exact  elasticity  theory, 

Eq.  (23),  are  also  given  in  this  table  for  zero-order  and  first-order  antisymmetric  Lamb 
waves,  respectively.  The  comparison  places  constraints  on  the  calculated  values  of  the 
correction  factors  at  low  and  high  frequencies. 


Table  Bl  —  Dispersion  Relation  for  Antisymmetric  Waves  in  the 
Limit  of  High  and  Low  Frequencies 


fed 

Smaller  Root 

Larger  Root 

Small 

72  =  K\y l 

Asymptote  for  kgd  at 

(exact  theory:  72  =  7 £) 

(fegd)J  =3  k]/k\ 

(exact  theory  at 
(fesd)J  =tt2/4) 

Large 

72  =  k2 

1  l 

-v? 

II 

M 

(exact  theory:  72  =  7^) 

(exact  theory:  7 2  =7^) 

Inserting  the  correction  factors  for  symmetric  waves  into  the  pertinent  equations  of 
motion,  Eqs.  (51)  and  (52),  leads  to  the  dispersion  relation 
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75  -  K-Wd  -H-lild  "  2> 

=  0.  (B6) 

Kl(yd  -2)  jk|(W)2  -j(kd)2 72  + 

The  smaller  root  of  this  quadratic  equation  is  identified  with  the  zero-order  symmetric 
Lamb  wave.  The  larger  root  can  be  identified  with  the  first-order  symmetric  Lamb 
wave.  There  is,  again,  a  vertical  asymptote  if  the  dispersion  curve  is  represented  in  the 
form  of  7  as  a  function  of  the  dimensionless  wavenumber  kgd.  Equation  (B6)  is  then 
transformed  into 

72  -  >cWd  -*l( y2d  -  2)y2 

=  0.  (B7) 

k-U yd-2)  -|^(fcsd)2  -±(ksd)2 72  +  K26ydy2 

By  solving  for  the  two  roots  of  this  equation,  one  can  again  derive  the  behavior  for 
low  and  high  frequencies  and  can  compare  the  results  with  the  exact  theory  for  sym¬ 
metric  waves.  This  is  shown  in  Table  B2.  If  one  computes  the  correction  factor  k2  ac¬ 
cording  to  Eq.  (Bl),  where  0  is  determined  by  0  =  (1/d)./^  (u/z)dz,  an  essential  difficulty 
arises.  Both  numerator  and  denominator  change  sign  at  intermediate  fed.  Although  the 
values  of  k23  for  low  and  high  frequencies  are  quite  acceptable,  the  sign  change  does  not 
occur  at  the  same  fed;  as  a  consequence,  the  correction  factor  goes  to  infinity  at  a  certain 
value  of  fed.  No  remedy  for  this  anomalous  behavior  has  been  found.  It  is  expected  that 
the  same  effect  may  occur  in  the  correction  factors  K4  and  k2  . 


Table  B2  —  Dispersion  Relation  for  Symmetric  Waves  in  the 
Limit  of  High  and  Low  Frequencies 


fed 

Smaller  Root 

Larger  Root 

Small 

2  __  2  2 

7  =  «47p 

Vertical  asymptote  at 

(exact  theory;  y2  =  72) 

(fesd)2  =  3K2by2d 
(exact  theory:  asymptote 
at  (fesd)2  =  7r2 ) 

Large 

72  =  *2 

72  =  K\y2d 

(exact  theory:  y2  =  7^) 

(exact  theory:  y2  =  7^) 
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QUADRATIC  APPROXIMATION  OF  DISPLACEMENT  COMPONENTS 


The  basic  mathematical  model  leading  to  thin-  and  thick-plate  theories  is  the  series 
expansion  of  the  displacement  components  u,  w  in  terms  of  the  coordinate  z,  perpendicular 
to  the  phases  of  the  plate.  So  far  no  higher  terms  than  linear  have  been  considered,  as  is 
shown  in  Eqs.  (27): 


u(x,y,z,t)  =  U(x,y,t)  +  z<t>x(x,y,t) 

v(x,y,z,t)  =  V(x,y,t)  +  z<py(x,y,t)  (27) 

w(x,y,z,t)  =  W(x,y,<)  +  zx(*,y,t)  . 


The  problem  encountered  with  the  correction  factor  in  symmetric  waves,  namely  that 
it  becomes  negative  for  low  fed,  prompted  an  effort  to  remedy  this  by  inclusion  of  higher 
order  terms.  Thus,  for  symmetric  waves  the  proposed  mathematical  model  for  the  dis¬ 
placement  coordinates  would  be 


u(x,z,t)  =  U(x,t)a+rz2/d2), 
w(x,z,t)  =  xz/d, 


(Cl) 


where  r  is  an  adjustable  constant.  The  shape  of  the  cross  section  with  the  assumptions  of 
Eq.  (Cl)  is  parabolic,  corresponding  to  the  visual  suggestion  of  the  symmetric  wave  de¬ 
picted  in  Fig.  Cl. 


ANTISYMMETRIC 


SYMMETRIC 


Fig.  Cl  —  Lamb  waves 
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Inserting  Eq.  (Cl)  into  the  equations  of  motion  of  Eq.  (37), 


t  bMxy 
dx  dy 


-Qx 


dMvx  dM  32LV 

_ yjL  + _ 3L-Q  =0 _ l 

dx  dy  y  dt2 


dQx  dQv  d2  Tz 

_ _  + _ y  -p — — , 

3jc  3y  dt2 


with  the  definitions  of  Eq.  (38) 


(37) 


leads  to  the  definition 


(C2) 


41 


DUBBKLDAY 


As  before,  there  is  some  latitude  in  the  choiee  of  the  new  term  rd//.  The  following  defini¬ 
tion  appears  plausible,  and  its  asymptotic  behavior  is  such  that  the  desired  value  for  k;  at 
high  frequency,  namely  7^,  as  obtained  before,  is  not  influenced: 


2  rdU 


(C3) 


This  definition  of  k\  is  an  improvement  in  the  sense  that  its  value  is  positive  over 
most  of  the  frequency  region.  The  same  problem  arises,  though,  as  was  encountered  with 
the  correction  factors  discussed  in  Appendix  B,  namely  a  sign  change  in  the  denominator 
that  causes  the  value  of  k*  to  go  to  infinity  at  a  certain  value  of  kd. 
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C  OMPARISON  WITH  OTHER  WORK  ON  APPROXIMATE  PLATE  THEORY 


In  the  literature  one  finds  other  work  on  the  subject  of  correction  factors  jn 
approximate  theories  for  thick  plates,  for  both  antisymmetric  and  symetric  waves.  It 
appears  that  the  distinction  between  such  papers  and  the  present  analysis  is  not  easily 
appreciated  [1)1].  The  present  study  follows  the  introduction  by  Mindlin  { D2 1  of  a  cor¬ 
rection  factor  for  effective  shear  modulus  and  proceeds  to  derive  an  analytic  expression 
for  this  correction  factor.  This  is  different  from  the  treatment  of  Mindlin,  who  determines 
the  correction  factor  by  fixing  the  high  frequency  limit  of  the  dis|>er.sion  relation  by  the 
Rayleigh  wave  speed.  In  this  study,  an  analogous  treatment  is  applied  to  symmetric 
waves  in  infinite  plates  and  a  corresponding  effective  shear  modulus  is  obtained  in  this 
case.  This  is  essentially  different  from  the  work  of  Kane  and  Mindlin  [ D 3 ]  and  Mindlin 
and  Medick,  [D4|  whose  analyses  apply  to  thickness  vibrations  in  finite  plates.  Some 
general  comments  on  the  classification  of  waves  in  plates  are  in  order,  to  clarify  the 
distinction. 

Solutions  of  partial  differential  equations,  including  the  wave  equation,  are 
determined  by  the  given  boundary  conditions.  A  natural  way  of  treating  waves  in  an  in¬ 
finite  solid  is  to  distinguish  dilatational  and  shear  waves,  since  this  reflects  the  general 
property  that  a  vector  field  can  be  represented  as  a  sum  of  an  irrotational  and  solenoidal 
field.  In  an  infinite  solid  the  two  waves  can  exist  independently.  In  the  study  of  waves 
propagating  parallel  to  the  faces  of  a  plate,  where  by  definition  the  thickness  of  the  plate 
is  small  compared  with  the  dimensions  parallel  to  the  plate,  the  dilatational  and  shear 
waves  are  coupled.  A  logical  division  in  this  case  is  that  of  antisymmetric  and  symmetric 
waves,  based  on  the  geometric  appearance  of  the  vibrating  plate.  The  general  theory  of 
Lamb  waves  is  developed  along  this  line  of  thought,  and  this  is  also  the  standpoint  of 
the  present  study  since  its  object  of  study  is  the  propagation  of  waves  in  extended  struc¬ 
tures.  Each  of  these  two  wave  types  is  composed  of  a  dilatational  part  and  a  shear  part. 

A  characteristic  phenomenon  in  Lamb  waves  is  the  appearance  of  higher  order  modes 
when  the  frequency  is  increased.  The  physical  explanation  of  these  higher  order  modes 
is  found  in  the  occurrence  of  standing  waves  in  the  direction  perpendicular  to  the  faces 
of  the  plate  that  are  coupled  to  the  traveling  waves  in  the  direction  parallel  to  the 
faces  of  the  plate.  These  thickness  modes  can  be  divided  into  thickness-stretch  and  thick¬ 
ness-shear  modes.  In  the  thickness-stretch  mode  the  particles  move  perpendicular  to  the 
faces  of  the  plate.  In  the  thickness-shear  mode  the  particles  move  parallel  to  the  faces 
of  the  plate.  Both  thickness  modes  may  occur  in  an  antisymmetric  and  a  symmetric 
version. 

One  could  obviously  start  one’s  analysis  by  emphasizing  the  thickness  modes  and 
treat  the  propagation  along  the  plate  as  a  secondary  effect.  This  viewpoint  is  represented 
in  a  review  of  plate  theory  by  Mindlin  [D5]  who  is  mostly  concerned  with  application 
of  the  general  theory  to  vibrations  of  crystal  plates.  It  is  the  description  of  choice 
in  the  analysis  of  vibrations  in  finite  plates,  in  contrast  to  the  study  of  propagation  of 
traveling  waves  in  infinite  plates.  Therefore,  it  is  clear  that  the  thickness  modes  of 
vibration  play  a  major  role  in  Refs.  D.i  and  D  l,  where  finite  plates  are  considered. 
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Oih*  may  obtain  an  overview  of  the  various  wave  types  in  plates  by  considering 
the  relationships  between  the  amplitudes  for  antisymmetric  ami  symmetric  waves, 
and  the  concomitant  dispersion  relations,  similar  to  the  discussion  in  Ref.  1)5.  The 
following  equations  are  identical  to  Kqs.  (22).  (23).  and  (21  I  in  the  mam  text 
I  he  equations  for  the  amplitudes  of  the  antisymmetric  waves,  under  homogeneous 
boundary  conditions,  are 


Bjk2+s2  )sinh(qd)  -  CQ  2 iks  sinh(sd)  =  0 
Ba  2 ikq  cosh (qd)  +  Ca(h2+s2  )  cosh (sd)  =  0, 
with  the  accompanying  dispersion  relation  for  antisymmetric  waves, 

(k2  +s2  )2  sinh(c/d)cosh(sd)  -  4k2  qs  cosh(qd  Isinhfsd  1  -  0 


(1)1  l 


(D2) 


The  equations  for  the  amplitudes  of  the  symmetric  waves,  under  homogeneous  boundary 
conditions,  are 


A  Ak2 +s2  )2  cosh(qd)  -  D  2iks  cosh(sd)  =  0 

*  S 

As  2 ikq  sinh(qd)  +  B>s(k2  +s2  )sinh(sd)  =  0, 
with  the  accompanying  dispersion  relation  for  symmetric  waves, 

(k2  +  s2  )2  cosh(<jd)sinh(sd)  -  4 k2  qs  sinh(qd)cosh(sd)  =  0. 


(D3) 


(D41 


the  lowest  frequency  at  which  thickness  modes  may  occur  is  determined  by  the 
condition  k  =  0.  This  means,  physically,  that  there  is  no  wave  propagation  in  the  direction 
of  the  faces  of  the  plate:  the  laminae  of  the  plate  move  in  unison,  in  the  direction  of 
the  faces  of  the  plate  for  thickness-shear  modes,  and  perpendicular  to  the  faces  for  thick¬ 
ness-stress  modes.  The  propagation  speeds  of  the  pertinent  waves  in  the  direction 
perpendicular  to  the  faces  are  the  shear  wave  speed  and  the  dilatational  wave  speed, 
respectively.  If  k  =  0.  it  follows  that  qd  =  ik^d  and  sd  =  /fe  d.  One  sees,  then,  that, 
given  I;  =  0,  the  following  combinations  are  solutions  to  Kqs.  (01)  through  (04). 
Antisymmetric  thickness-shear  is  obtained  for  B  -  0,  cos (A*?di  =  0,  or  w  =  (2n  +  1 ) 

( -  2)cs/d.  Antisymmetric  thickness-stretch  modes  are  obtained  for  CQ  =  0,  sin(fcrfd)  =  0. 
or  ca  =  niTcd/d.  Symmetric  thickness  shear  modes  occur  for  A  =  0,  sini/o.d)  =  0,  or 
w  =  nrc$ld.  Symmetric  thickness-stretch  modes  occur  for  D  '=  0.  cos(/rjd)  =  0,  or 
sc  =  (2n+l  )(t:/2  )('d/d.  Here  n  =  1,2,  3...;  cs  is  the  shear  wave  speed;  and  cd  is  the 
dilatationa.  wave  speed. 


The  structure  of  the  Kqs.  (01)  through  I  04)  leads  to  the  question  of  what  type  of 

waves  will  result  if  the  condition  k  -■  0  is  replaced  by  the  condition  I;2  +  s2  =  0. 

Since  s2  =  k2  -  k2.  it  follows  that  r2  =  2 c2  for  this  type  of  wave.  Given  this  condition, 

the  first  possibility  for  antisymmetric  waves  is  that  B  -  0,  and  sin(fr  d/  /2)  =  0,  or 

nit  T>c  ,  d  This  type  of  wave  is  indicated  by  the  term  equuoluminnl  by  Mindlin  1 01  ] , 

■'ince  the  amplitude  of  the  dilatational  part  of  the  wave  is  zero  in  this  case.  Actually  the 

thickness-shear  mode  des'-nbed  before  would  equally  well  qualify  for  this  name  for  the 

■amr  reason.  The  analogous  complementing  case  can  be  found  by  setting  C  =  0  and 

eodiiuTi  0.  The  latter  condition  cannot  be  satisfied,  since  here  a2  0.5/r  -  /,-2,  which 

X  (I ' 
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is  always  positive.  Therefore  this  case  has  no  physical  reality.  In  a  similar  way  a  symmet¬ 
ric  “equivoluminal”  wave  may  be  derive  by  setting  Ag  =  0,  with  cosh(sd)  =  o,  which 
leads  to  oj  =  (2o+l  )(7r/2)  J~2cs/d. 

In  Ref.  D5  Mindlin  proposes  a  correction  factor  k  to  improve  the  correspondence 
between  approximate  plate  theory  and  exact  theory.  Mindlin  develops  his  theory  using 
the  classical  power-series  expansion  method  that  originated  with  Poisson  and  Cauchy. 

He  adjusts  the  constant  k  in  such  a  way  that  the  first  appearance  of  the  thickness  modes 
occurs  at  the  correct  frequency.  This  leads  to  the  value  of  rr2  / 1 2  for  his  correction  fac¬ 
tor,  for  both  antisymmetric  and  symmetric  waves.  Since  the  displacement  components 
in  thickness  modes  vary  in  a  harmonic  way  as  a  function  of  the  coordinate  perpendicular 
to  the  plate,  it  would  follow  that  expansion  in  terms  of  a  Fourier  series  rather  than  a 
power  series  promises  a  better  correspondence  between  the  approximate  and  exact  theo¬ 
ries.  Employing  a  Fourier  series  expansion,  however,  is  found  to  result  in  poor  behavior 
of  waves  traveling  parallel  to  the  face  of  the  plate  at  high  frequencies.  Again,  the  latter 
circumstance  is  of  little  importance  if  the  object  of  study  is  vibrations  of  finite  plates, 
as  in  Refs.  D3  and  D4. 

In  an  earlier  publication  [D2]  Mindlin  had  introduced  a  similar  correction  factor  to 
improve  the  high-frequency  behavior  of  the  dispersion  relation  of  the  zero-order  antisym¬ 
metric  mode  (flexural  wave)  in  approximate  plate  theory.  There,  he  so  fixes  the  factor  k 
that  at  high  frequency  the  wave  speed  approaches  the  Rayleigh  wave  speed.  Mindlin  com¬ 
ments  that  thus  a  compromise  has  to  be  found  between  this  high-frequency  adjustment  of 
k  and  the  low-frequency  adjustment,  which  amounted  to  setting  k  equal  to  u2/\/l2,  as  dis¬ 
cussed  before. 

Here  an  important  difference  appears  between  the  present  study  and  the  referenced 
papers.  In  the  present  work,  the  emphasis  is  on  wave  propagation  along  an  infinite  plate, 
represented  by  zero-order  antisymmetric  and  symmetric  waves,  which  correspond  to  the 
smaller  root  of  the  quadratic  dispersion  relation  for  flexural  and  extensional  waves  in 
thick-plate  theory.  Moreover  the  correction  factor  in  the  present  report  is  not  fixed  at 
either  the  high-frequency  or  low-frequency  limit,  but  is  found  by  comparison  of  the 
integrals  in  the  approximate  theory  with  those  obtained  from  exact  elasticity  theory  over 
the  whole  frequency  range.  No  “compromise"  is  needed,  because  here  k  is  frequency 
dependent,  and  moreover  the  low-frequency  behavior  of  the  zero-order  mode  (the  smaller 
root  of  the  quadratic  equation)  is  independent  of  the  value  of  k.  Thus  both  the  high- 
frequency  and  low-frequency  behaviors  are  properly  accounted  for  in  the  present  analysis. 

Kane  and  Mindlin  [D3]  analyze  the  behavior  of  extensional  vibrations  in  plates  at 
high  frequencies.  One  might  surmise  that  their  work  is  analogous  to  Mindlin’s  study  of 
flexural  motions  in  plates  [D2]  in  the  sense  that  a  comparable  treatment  is  applied  to 
symmetric  waves,  but  that  is  not  the  case.  The  study  of  Kane  and  Mindlin  is  limited  to 
vibrations  of  finite  plates,  and  as  a  consequence  they  fix  the_cqirection  factor  k  at  the 
low-frequency  end  by  assigning  to  it  the  familiar  value  n2 1\/ 12.  An  unfortunate  aspect 
of  the  basic  approach  of  these  authors  is  that  the  pertinent  correction  factor  is  intro¬ 
duced  into  the  differential  equations  of  motion  of  elasticity  theory  (Eq.  (20)  in  Ref. 

D3).  There  is  no  valid  reason  to  modify  the  fundamental  equations  of  elasticity  theory 
used  all  through  the  literature.  Instead,  the  correction  factor(s)  should  properly  appear 
in  the  plate-stress  equations,  where  approximations  are  introduced  for  the  variation  of 
the  displacement  components  as  a  function  of  the  coordinate  perpendicular  to  the  plate. 
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Mmdlin  and  Mediek  (1)1 1  expand  on  the  ideas  of  Kane  and  Mindlm  m  several 
respects.  In  the  first  plare,  instead  of  using  a  straight  power-series  expansion,  they  use 
an  expansion  in  Legendre  polynomials  for  the  displacement  components.  They  expect 
that  such  a  series  may  improve  the  quality  of  the  approximation  and  may  simplify  the 
task  of  assessing  the  effect  of  truncation.  In  the  second  place,  different  correction  fac¬ 
tors  are  introduced  that  are  similar  to  the  correction  factors  of  the  present  study.  How¬ 
ever,  these  authors  also  make  the  point  that  the  correction  factors  can  be  fixed  only  at 
a  finite  number  of  points  in  the  frequency  range  in  order  to  improve  the  correspondence 
between  approximate  and  exact  theories  with  respect  to  the  dispersion  curves.  Quite 
emphatically,  Mmdlin  and  Mediek  opt  for  improving  the  thick-plate  theory  at  the  limit 
of  zero  kd,  in  order  to  improve  the  treatment  of  finite  plates,  for  which  the  thickness 
modes  are  the  most  important. 

In  connection  with  the  above  discussion,  the  following  point  should  be  stressed. 

In  the  study  of  waves  in  plates,  one  may  concentrate  on  the  propagation  of  traveling 
waves  along  the  faces  of  infinite  plates  or  on  the  standing  waves  perpendicular  to  the 
faces  of  the  plates  for  finite  plates.  The  two  cases  correspond  to  the  smaller  and  larger 
roots  of  the  quadratic  equations  expressing  the  dispersion  relations  in  thick-plate  theory. 
Considering  the  fact  that  the  character  of  the  motion  in  the  two  cases  is  quite  different, 
there  is  no  reason  to  believe  that  one  and  the  same  correction  factor  would  suffice  in 
both  cases  over  the  whole  frequency  range.  Several  correction  factors  are  needed.  Their 
role  in  connection  with  the  low-  and  high-frequency  behaviors  is  listed  in  Tables  D1  and 
D2.  One  sees  in  Table  D1  that  the  correc  tion  factor  for  the  shear  modulus  in  antisym¬ 
metric  waves  X!  does  not  influence  the  smaller  root  at  low  frequency.  If  one  assumes 
that  the  factor  x for  the  larger  root  has  a  limit  of  1  at  the  low-frequency  limit,  the  low- 
frequency  limit  of  x,  for  this  branch  will  be  the  familiar  -2\/12.  Remember,  though, 
that  this  factor  is  computed  here  on  the  basis  of  that  branch  of  the  dispersion  relation 
that  corresponds  to  first-order  antisymmetric  Lamb  waves.  Table  D2  shows  that  for 
symmetric  waves  the  correction  factor  governing  the  behavior  of  thickness  vibrations 
at  low  frequency  (x6  )  is  not  even  the  same  nominal  factor  as  the-  one  that  determines 
the  high-frequency  limit  for  traveling  waves  (x2  in  the  smaller  root). 


Table  D1  —  Dispersion  Relation  for  Antisymmetric  Waves 
in  the  Limit  of  High  and  Low  Frequencies 


Smaller  Root 

Larger  Root 

Small  kd 

7J  =  k)  y\ 

Asymptote  for  k^d  at 

(exact  theory:  y7  =  y b2  ) 

=  .V;/k' 

(exact  theory  at 

(Ay l)7  -  'll 

Large  kd 

y2  =  k] 

(exact  theory:  y7  =  y]{) 

(exact  theory:  y  2  y  ],  1 
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Table  D2  —  Dispersion  Relation  for  Symmetric  Waves  in  the 
Limit  of  High  and  Low  Frequencies 


Smaller  Root 

Larger  Root 

Small  kd 

2  _  2  2 
y*  r  7 

4  '  p 

Vertical  asymptote  at 

(ksd)7  =  3  k7  y7d 

(exact  theory:  y7  =  7*  ) 

(exact  theory:  asymptote 

at  (kgd)7  =  7t 7 ) 

Large  kd 

r«  n 

* 

II 

CN 

?** 

2  2  2 

7  yd 

(exact  theory:  y7  =  y7R  ) 

(exact  theory:  y7  =  y7R  ) 

In  conclusion,  the  present  report  is  distinguished  from  Refs.  D2  through  D4  in  the 
following  respects. 

1.  It  presents  a  thick-plate  theory  for  plates,  for  which  the  dimensions  parallel  to 
the  plate  are  considerably  larger  than  the  thickness,  at  such  frequencies  that  the  zero- 
order  antisymmetric  and  symmetric  Lamb  waves  are  the  only  ones  of  importance. 

2.  It  gives  a  method  of  calculating  the  correction  factor  for  effective  shear  modulus 
in  antisymmetric  waves  by  comparing  thick-plate  theory  with  the  results  of  exact  elasti¬ 
city  theory,  instead  of  fixing  the  correction  factor  at  one  point  in  the  frequency  range. 

3.  It  presents  a  thick-plate  theory  for  symmetric  waves  with  a  correction  factor  for 
effective  shear  modulus,  corresponding  to  thick-plate  theory  of  antisymmetric  waves;  it 
also  gives  a  method  of  computing  this  factor  by  comparing  thick-plate  theory  with  the 
results  of  exact  elasticity  theory,  instead  of  fixing  the  factor  at  one  point  in  the  frequency 
range. 
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